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PART  I  -  CAVITATION  INCEPTION 


X 


I.  INTRODUCTION 


Background 

As  a  consequence  of  the  equation  of  state  one  of  the  important 
differences  between  dynamics  of  liquids  and  gases  is  that,  liquids  can 
withstand  tensile  stresses  whereas  gases  cannot.  If  sufficient  tensile 
stresses  are  applied,  the  liquid  ruptures  to  form  cavities.  Theoretical 
predictions  indicate  that  pure  homogeneous  liquids  can  withstand  con¬ 
siderable  amount  of  tensile  stresses  of  the  order  of  thousands  of 
atmospheres^.  So  far  carefully  treated  liquid  samples  are  known  to 
withstand  only  moderate  tensile  stresses  of  the  order  of  hundreds  of 
atmospheres^1).  The  reason  for  this  discrepancy  has  led  to  the  postu¬ 
lation  of  presence  of  impurities  in  the  liquid  samples  commonly  known 
as  "nuclei".  The  initiation  of  vapor  growth  on  these  nuclei  under  the 
action  of  the  tensile  stress  to  result  in  visible  cavity  is  known  as  the 
inception  of  cavitation.  It  is  clear  that  if  these  cavities  once  formed 
are  subjected  to  compressive  stress,  they  would  tend  to  disappear.  In 
general,  cavitation  is  the  study  of  appearance,  development  and  dis¬ 
appearance  of  cavities  in  a  body  of  liquid. 

If  the  cavities  contain  mainly  the  vapor  form  of  the  cavitating 
liquid,  the  process  is  known  as  "vaporous"  cavitation;  on  the  other 
hand  if  they  mainly  contain  dissolved  gas  which  is  one  of  the  most 
common  forms  of  impurity  in  a  liquid,  it  is  termed  as  "gaseous"  cavi¬ 
tation.  Cavitation  is  of  common  occurrence  on  ship  propellers, 
hydraulic  machinery  such  as  liquid  primping  devices  and  on  hydrofoils 
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used  as  lifting  devices  for  hydrofoil  craft.  In  almost  all  cases  occur¬ 
rence  of  cavitation  is  an  undesirable  phenomenon  due  to  limitations  on 
the  design  of  liquid  flow  devices,  the  potential  of  damage  to  the  solid 
surfaces  of  the  body,  unsteady  forces  and  noise. 

Since  the  onset  of  cavitation  distinguishes  the  limit  for 
cavitation-free  operation,  it  has  been  the  subject  of  considerable 
experimental  research.  In  this  type  of  work  it  is  of  great  importance 
to  predict  the  onset  or  inception  of  cavitation  on  prototypes  from  model 
studies.  With  regards  to  this  prediction,  scaling  laws  which  control 
the  behavior  of  flows  under  cavitating  conditions  have  to  be  found.  There¬ 
fore  by  "scaling"  it  is  meant  to  extrapolate  the  cavitation  behavior 
from  one  set  of  conditions  to  another.  To  achieve  the  goal  of  finding 
appropriate  scaling  laws,  considerable  experimental  research  has  been 
conducted  on  certain  simple  bodies.  With  these  studies  it  is  hoped  that 
a  clear  understanding  of  the  basic  mechanism  of  cavitation  is  obtained 
in  order  to  predict  the  behavior  of  flows  under  cavitating  conditions  in 
complex  situations. 

One  of  the  important  parameters  in  describing  cavitating  flow  is 
the  cavitation  number  c  defined  as 

o  -  (Pq-Pv  )/  ■£  P  Uq  where, 

Pq  is  the  reference  pressure,  Uq  is  the  reference  velocity,  p  is  the 
liquid  density  and  P  is  the  vapor  pressure  of  the  liquid  at  its  bulk 
temperature.  The  value  of  o  at  which  inception  of  cavitation  occurs  is 
designated  as  o^.  In  water  tunnel  experiments  Pq  and  Uq  refer  to  the 


tunnel  static  pressure  and  velocity  respectively,  measured  some 
distance  away  from  the  model.  In  most  water  tunnels,  experimental 
conditions  of  cavitation  inception  are  achieved  by  lowering  Pq  until 
signs  of  cavitation  are  visible  in  the  neighborhood  of  the  test  model 
while  the  tunnel  velocity  is  held  fixed.  If  a  is  lowered  sufficiently 
below  that  of  o^  a  developed  cavity  is  formed  on  the  model  surface  and 
if  the  tunnel  static  pressure  is  increased  slowly,  cavitation  disappears 
and  condition  of  existence  of  last  traces  of  cavitation  on  the  model  is 
known  as  "desinent"  cavitation.  The  corresponding  a  is  designated  by 
Oj  and  has  been  used  to  describe  the  limit  of  the  noncavitating  flow 
regime  notably  by  Holl^.  In  general  it  has  been  found  that  Oj  >  Oj  and 
this  can  be  viewed  as  a  hysteresis  phenomenon.  But  o  is  not  the  only 
pertinent  physical  parameter  in  determining  cavitation  inception,  and 
a  complete  set  of  parameters  determining  in  general  is  yet  to  be 
found^*  An  example  of  this  is  illustrated  in  measurements^  of 
cavitation  inception  index  on  different  sized  cylindrical  bodies  with 
hemispherical  noses.  Some  of  the  results  from  Ref.  (4)  are  shown  in 
Fig.  1.  Oj  is  seen  to  vary  with  both  size  of  the  body  and  tunnel  velocity. 
This  dependence  of  Oj  on  the  size  of  the  body,  velocity  and  other 
variables  in  determining  inception  of  cavitation  is  commonly  termed 
the  "scale  effect".  Additional  inception  tests  conducted  on  different 
types  of  streamlined  bodies  such  as  a  1.5  calibre  ogive^  and  various 
hydrofoil  sections^’  ^  have  exhibited  different  effects  with  change  of 
scale. 

Cavitation  inception  studies  are  normally  carried  out  in  the 
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many  different  tunnels  throughout  the  world.  Eich  of  these  possesses 
different  hydrodynamic  and  physical  characteristics.  A  program  to 
carry  out  a  comparison  of  cavitation  inception  measurements  in  many 
of  these  tunnels  was  recently  organized  by  the  Cavitation  Committee  of 
the  International  Towing  Tank  Conference  (I.  T.  T.  C. ).  The  purpose  of 
the  investigation  was  to  determine  the  incipient  cavitation  index  a j  on  a 
particular  axi-symmetric  half  body  (see  Fig.  In  this  program  each 
laboratory  was  advised  to  use  its  own  normal  routine  of  performing 
inception  measurements;  no  suggestions  were  given  regarding  the 
diameters  of  the  model,  the  way  of  measuring  tunnel  static  pressure, 
water  speed  and  so  on.  Other  experimental  details  may  be  found  in 
Ref.  (7).  Some  of  the  results  from  this  study  are  shown  in  Fig.  2  from 
which  it  is  evident  that  considerable  differences  in  the  values  of  ob¬ 
served  a.  were  found  in  different  tunnels.  It  is  also  discussed  in  Ref. 
(7)  that  the  appearance  or  form  of  the  cavitation  could  be  different  in 
the  various  facilities.  One  of  the  principle  design  differences  between 
various  tunnel  facilities  is  the  presence  or  absence  of  a  "resorber" 
in  the  tunnel  circuit.  The  purpose  of  a  resorber  is  to  put  into  solution 
the  free  air  bubbles  created  during  the  process  of  tunnel  operation. 

This  is  achieved  by  subjecting  the  flow  to  high  pressures  before  re¬ 
entering  the  test  section.  Tunnel  facilities  which  do  not  contain  a 
resorber  tend  to  have  an  abundance  of  macroscopic  air  bubbles  in  the 
water  approaching  the  test  body  (readily  visible  under  stroboscopic 
illumination).  These  macroscopic  bubbles  when  passing  through  the 
low  pressure  region  in  the  neighborhood  of  the  model  explode  to 
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vaporous  cavities  and  subsequently  collapse.  Photographs  of  this  type 
of  cavitation  have  been  obtained  by  Knapp'  '.  A  different  type  of 
inception  is  observed  in  tunnel  facilities  which  do  have  a  resorber. 

This  type  has  been  described  by  Acosta^,  "At  this  stage  we  would 
like  to  describe  the  cavitation  that  appears  on  the  I.  T.  T.  C.  test  body 
in  the  high  speed  water  tunnel.  In  practically  all  of  the  present  tests 
no  free  bubbles,  or  very  few,  in  the  working  section  upstream  of  the 
test  body  were  observed  even  up  to  speeds  of  20.0  meters  per  second. 
This  is  in  cortrast  to  many  other  water  tunnels  in  which  a  copious  num¬ 
ber  of  free  air  bubbles  are  plainly  evident  —  for  example,  the  Naval 
Ordnance  Test  Station  Vertical  Water  Tunnel.  Furthermore,  as  the 
point  of  cavitation  inception  on  the  body  is  reached,  no  minute  traces 
of  cavitation  bubbles  were  seen  on  the  I.  T,  T.  C.  test  body.  In 
practically  every  case  when  cavitation  inception  occurred,  it  took  place 
suddenly  and  abruptly  in  the  form  of  a  small  band  all  around  the  nose  of 
the  body.  This  appearance  was  characteristically  abrupt  rather  than 
gradual  and  there  were  no  "precursors"  suggesting  that  cavitation 
inception  was  imminent' .  As  further  discussed  in  Ref.  (9)  this  form  of 
cavitation  is  termed  "band"  inception.  But,  on  the  same  I.  T.  T.  C.  test 
body  at  higher  speeds  another  type  of  cavitation  inception  was  also 
observed  and  is  termed  as  "spot  cavitation",  which  commenced  at  a 
pinpoint  or  spot  and  subsequently  spread  out  in  the  form  of  a  "V"  rather 
like  a  wake.  A  photograph  of  "band"  type  inception  as  observed  in  Ref. 
(9)  is  shown  in  the  bottom  photograph  of  Fig.  1  of  Part  II. 

Now  called  the  Naval  Undersea  Research  and  Development  Center. 
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Some  of  the  differences  observed  in  the  values  of  a.  in  the 

different  facilities  shown  in  Fig.  2  may  be  attributed  to  the  different 

types  of  cavitation  inception  just  discussed.  Indeed,  the  work  done  at 

D.  T.  M.  B.  shows  that  values  of  o.  measured  in  a  36  inch  tunnel  with  a 

i 

resorber  were  systematically  lower  by  about  15-20  percent  than  the 
values  of  measured  in  the  same  tunnel  without  the  resorber  feature. 
The  I.  T.  T.  C.  report  '  '  was  not  able  to  resolve  these  differences  as 
well  as  others  on  geometrically  similar  bodies  tested  in  different 
laboratories.  One  may  conclude  from  this  work  that  the  test  body  used 
or  procedures  adopted  did  not  lead  to  a  useful  means  of  predicting 
cavitation  inception,  or  to  the  ability  of  translating  the  results  of  one 
experimental  facility  to  another. 

Some  Physical  Features  of  Cavitation  Inception:  It  has  been  mentioned 
before  that  "nuclei"  or  weak  spots  are  necessary  to  initiate  the  phase 
change  that  results  in  the  cavitation  that  is  observed  at  inception.  It  is 
apparent  that  in  considering  flows  about  a  submerged  body  as  in  the 
experiments  mentioned  before,  the  nuclei  can  come  from  two  sources; 
namely,  from  within  the  free  stream,  and  the  surface  of  the  test  body 
itself.  Nuclei  of  the  first  type  are  termed  as  "stream  nuclei"  and  the 
latter  as  "surface  nuclei".  The  physical  nature  of  the  free  stream 
nuclei  have  been  hypothesized  to  be  of  the  form  of  free  gas  bubble 8^^’  ^ 
gas  trapped  in  the  crevices  of  solid  particles  in  the  free  stream^  or 
solid  particles  themselves^  ^  which  are  present  in  the  fluid.  It  has 
Now  the  Naval  Ship  Research  and  Development  Center. 
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Some  of  the  differences  observed  in  the  values  of  a.  in  the 
different  facilities  shown  in  Fig.  2  may  be  attributed  to  the  different 
types  of  cavitation  inception  just  discussed.  Indeed,  the  work  done  at 
D.  T.  M.  B.  shows  that  values  of  a.  measured  in  a  36  inch  tunnel  with  a 
resorber  were  systematically  lower  by  about  15-20  percent  than  the 

values  of  measured  in  the  same  tunnel  without  the  resorber  feature. 
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been  suggested'  ’  that  surface  nuclei  consist  of  gas  trapped  in  micro¬ 
scopic  crevices  of  the  body  surface.  To  understand  the  role  of  nuclei 
in  cavitation  inception  it  is  helpful  to  consider  the  static  stability  of  a 
stationary  spherical  gas  bubble  containing  the  vapor  of  the  fluid  and  a 
small  amount  of  permanent  gas.  It  is  found  from  the  results  of  such  an 
analysis  (see  for  example,  Strasberg^^)  that  for  a  given  radius  Rq 
there  is  an  ambient  pressure  below  which  the  bubble  becomes  unstable. 
This  pressure  is  designated  as  the  "critical"  pressure  Pc  which  depends 
on  tiie  initial  radius  Rq.  When  the  ambient  pressure  is  less  than  Pc,  the 
bubble  may  grow  under  the  condition  of  constant  vapor  pressure  and  gas 
content.  In  this  event  the  bubble  may  be  said  to  grow  by  vaporous 

cavitation  (i.  e, ,  gaseous  diffusion  does  not  occur).  It  can  be  shown  that 

-2 

as  Rq  becomes  larger  than  about  10  inches,  the  critical  pressure  for 
the  bubble  growth  (or  vaporous  cavitation)  approaches,  but  never  ex¬ 
ceeds,  the  vapor  pressure  i.e. ,  the  critical  pressure  is  always  less 
than  the  vapor  pressure.  Smaller  bubbles  with  radii  of  the  order  of 
2  X  10  inches  for  example,  would  have  a  critical  pressure  approxima¬ 
tely  one  atmosphere  less  than  the  vapor  pressure.  A  bubble  of  this 
size  in  water  at  room  temperature  would  therefore  require  a  tension  of 
this  magnitude  for  instability  to  occur.  Thus  we  can  see  that  for  any 
given  liquid  pressure  lower  than  the  vapor  pressure,  a  bubble  of  radius 
less  than  a  certain  minimum  radius  will  not  grow  or  cavitate  vaporously. 

This  minimum  radius  may  be  designated  as  the  "critical  radius"  R  . 

c 

We  would  therefore  expect  from  these  examples  that  if  a  sufficient 
number  of  nuclei  are  present  in  a  flowing  stream  with  radii  larger  than 


about  1(j  inches,  vaporous  cavitation  would  begin  as  soon  as  the 
pressure  in  the  neighborhood  of  a  submerged  body  falls  to  the  vapor 
pressure.  In  any  case,  since  the  minimum  pressure  on  the  body  is  the 
lowest  pressure  available  in  the  flow  it  has  been  customary  to  predict 
o\  by  equating  it  with  the  minimum  pressure  coefficient  of  the  body, 

-  cpmjn*  But  from  comparative  studies  on  the  I.  Tt  T.  C.  test  body  it 

may  be  observed  that  Cj  can  be  considerably  different  from  -  cpmjn  as 

can  be  seen  in  Fig.  2.  For  the  cases  a,<  -c  .  it  may  be  concluded 

ip  mm 

either  that  the  value  of  -  Cpmjn  was  not  known  correctly  or  that  the 
critical  pressure  of  liquid  was  not  vapor  pressure.  It  can  also  be  seen 
from  Fig.  2  that  some  of  the  values  of  o j  are,  in  fact,  greater  than 

-  Cp min>  indicating  occurrence  of  cavitation  at  pressures  higher  than 
vapor  pressure.  This  observation  suggests  that  a  process  other  than 
vaporous  cavitation  is  responsible  and  indeed  this  is  attributed  to 
gaseous  diffusion  and  is  discussed  by  Holl^  in  his  review  article. 

This  process  is  termed  "gaseous  cavitation"  in  contrast  to  vaporous 
cavitation.  Gaseous  cavitation  can  occur  at  any  ambient  pressure 
(unlike  vaporous  cavitation)  provided  that  the  liquid  is  supersaturated 


in  respect  to  the  partial  pressure  of  the  gas  within  the  bubble.  The 

dynamics  of  the  vaporous  and  gaseous  growth  of  a  nucleus  is  of  different 

time  scales.  Typical  times  for  vaporous  growth  are  of  tl*  order  of 

milliseconds'  whereas  for  gaseous  growth  they  may  be  of  the  order 
(15) 

of  few  seconds'  \  This  would  seem  to  preclude  gaseous  cavitation  as 


being  technically  important  because,  if  bubbles  move  with  the  free 
stream  they  do  not  remain  in  low-pressure  regions  long  enough  for 
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appreciable  growth.  However,  if  the  bubbles  are  trapped  on  a  fixed 
boundary,  sufficient  time  may  be  available  for  diffusion.  It  may  be 
imagined  that  both  processes  occur  in  cavitation  on  real  bodies,  as  will, 
in  fact,  be  seen  in  the  experiments  to  be  discussed  subsequently. 

Thus  cavitation  inception  process  can  be  assumed  to  be  depen¬ 
dent  on  two  main  factors,  namely,  (a)  sources  of  cavitation  nuclei,  and 
(b)  the  growth  potential  of  the  nucleus.  Following  are  some  detailed 
considerations  of  the  above  mentioned  two  factors. 

Sources  of  Cavitation  Nuclei:  As  mentioned  earlier  depending  on  their 
origin;  nuclei  can  be  classified  into  either  surface  nuclei  or  stream 
nuclei.  In  recent  years,  there  have  been  several  investigations  which 
have  shown  that  inception  of  cavitation  is  influenced  by  surface  charac¬ 
teristics,  which  presumably  control  the  characteristics  of  surface 
nuclei.  Holl  and  his  co-workers^’  have  carried  out  cavitation 

inception  studies  on  hemispherical  nosed  bodies  having  surfaces  of 
smooth  stainless  steel,  sandblasted  steel,  waxed  stainless  steel,  teflon, 
nylon  and  glass.  It  was  observed  that  teflon  was  easy  to  cavitate 
whereas,  glass  was  hard  to  cavitate.  Such  observations  prompted 
Holl'  ’  to  suggest  that  porosity  might  be  an  important  factor  influencing 
the  characteristics  of  surface  nuclei.  It  is  not  clear  from  their  findings 
whether  the  change  in  inception  characteristics  was  due  to  the  above 
mentioned  reason  or  due  to  the  change  in  the  nature  of  the  boundary  layer 
growth  on  bodies  having  different  surface  characteristics.  Peterson^ ^ 
recently  has  found  that  if  the  test  body  is  cleaned  with  acetone,  then 
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subjected  to  a  pressure  of  19,  000  psi  (presumably  effect  of  this  is  to 
drive  some  of  the  air  trapped  in  the  model  crevices  into  solution,  thus 
changing  effective  size  of  the  surface  nuclei)  this  produced  sixteen  per¬ 
cent  reduction  in  the  observed  value  of  a.  compared  with  the  values 
obtained  without  any  surface  treatment.  Another  interesting  experiment 
concerning  surface  nuclei  was  carried  out  by  Acosta  and  Hamaguchi'  \ 

In  their  experiments  the  nose  of  the  test  body  (the  same  test  body  used 
in  the  I.  T.  T.  C.  comparative  study)  was  dipped  into  a  mixture  of  a 
silicon  fluid  and  solvent.  The  object  of  this  procedure  was  to  coat  the 
surface  of  the  body  with  a  fluid  which  had  the  ability  to  dissolve  a  large 
amount  of  air  as  silicone  fluids  are  known  to  do.  The  inception  studies 
showed  that  when  the  silicone  coating  was  fresh  and  presumably  saturated 
with  air  incipient  cavitation  numbers  to  be  systematically  higher  than 
that  obtained  for  the  case  of  pure  water.  But  as  the  coating  got  old  and 
deaerated  due  to  evolution  of  cavitation  bubbles  the  values  of  o.  were 
found  to  be  lower  than  that  observed  for  the  case  of  pure  water.  It 
should  be  mentioned  that  the  tests  were  preliminary  in  nature,  but  they 
do  indicate  the  possibility  of  existence  of  surface  nuclei  in  the  form  of 
air  trapped  in  the  model  crevices.  Still,  the  knowledge  regarding  the 
role  of  surface  nuclei  in  the  mechanism  of  cavitation  inception  is 
obscure. 

The  role  of  free  stream  nuclei  in  the  form  of  free  air  bubbles 
in  the  process  of  cavitation  inception  has  been  studied  by  Ripken, 

Killen  and  Schiebe^®'  It  is  expected  from  the  stability  analysis 
of  gas  bubbles  considered  earlier,  that  the  si zo  spectrum  of  gas  bubbles 
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present  in  the  stream  would  be  of  importance  in  determining  the  value 
of  o..  An  acoustic  method  of  obtaining  such  a  spectrum  for  recirculat¬ 
ing  facilities  has  been  successfully  used  by  Schiebe^^.  Major  findings 
of  Ripken  and  Killen  summarized  in  Ref.  (20)  are  that  the  free  gas 
content  (in  the  form  of  gas  bubbles)  is  at  least  a  function  of:  a)  total 
gas  content,  which  is  the  sum  of  the  free  gas  content  and  the  dissolved 
gas  content,  b)  tunnel  pressure,  and  c)  tunnel  velocity.  They  also 
mention  that  the  water  tunnel  pump  is  probably  the  main  source  of  gas 
bubbles  and  that  the  free  gas  spectrum  would  be  a  unique  property  of 
individual  tunnel  facilities.  It  is  clear  from  earlier  considerations  that 
measurement  of  free  gas  spectrum  would  be  of  greater  importance  for 
tunnel  facilities  without  the  presence  of  a  resorber  than  for  facilities 
with  the  presence  of  a  resorber. 

Growth  Potential  of  the  Nucleus:  The  ability  of  the  available  nucleus  to 
grow  to  visible  size  would  definitely  be  influenced  by  the  magnitude  and 
duration  of  the  low  pressure  experienced  by  the  nucleus.  Photographic 
investigations  conducted  by  Knapp  and  Hollander '  '  years  ago  have 
indicated  that  the  growth  history  of  cavitation  bubbles  orginating  from 
free  stream  nuclei  is  only  slightly  influenced  by  the  presence  of  the 
boundary  layer  on  the  model  surface.  The  development  of  surface 
nuclei  into  visible  cavitation  would  be  expected  to  be  greatly  influenced 
by  the  presence  of  the  boundary  layer  on  the  model  surface  contrary  to 
the  behavior  of  free  stream  bubble  cavitation.  It  is  this  point  which  has 
received  very  little  attention  in  the  understanding  of  cavitation  inception 


mechanism.  During  Kermeen's'  '  experiments  on  hemispherical  nosed 
bodies,  it  was  observed  that  cavitation  inception  took  place  in  the  form 
of  a  ring  of  macroscopic  bubbles  around  the  nose  of  the  headform  and 
photographs  indicated  the  presence  of  microscopic  bubbles  with  dia¬ 
meters  smaller  than  typical  boundary  layer  thicknesses  ahead  of  the 
macroscopic  ones.  These  observations  prompted  Parkin  and  Kermeen1  ' 
to  study  the  influence  of  the  presence  of  the  boundary  layer  on  the 
mechanism  of  cavitation  inception.  Thus  from  their  now  classic  experi- 
ments'  '  of  1953  they  concluded  that  the  presence  of  a  boundary  layer 
plays  a  significant  role  in  the  process  of  cavitation  inception  on  stream¬ 
lined  bodies.  Their  results  can  be  summarized  as  follows.  The 
macroscopic  cavitation  visible  to  the  unaided  eye  was  sustained  by  very 
small  '•microscopic"  bubbles  which  grew  in  the  boundary  layer.  The 
microscopic  bubbles  grew  for  a  time  at  a  fixed  position  on  the  body 
upstream  from  the  macroscopic  cavitation.  These  bubbles  after  attain¬ 
ing  a  radius  of  approximately  one  half  the  displacement  thickness  of  the 
boundary  layer  were  then  entrained  into  the  flow  and  travelled  down¬ 
stream  to  be  distorted  into  irregular  patches  which  sustained  the 
macroscopic  cavitation.  After  further  reduction  of  tunnel  static 
pressure  the  evolution  of  microscopic  cavitation  bubbles  increased 
until  a  clear  attached  glassy  cavity  was  formed  on  the  model  surface. 
When  this  occurred  the  microscopic  cavitation  just  described  was 
completely  suppressed.  It  was  also  concluded  that  latter  stages  of  the 
growth  of  microscopic  bubbles  was  significantly  influenced  by  gaseous 
diffusion.  It  should  be  mentioned  that  flow  visualization  of  the  boundary 
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layer  was  done  by  introducing  a  dye  streak.  In  their  final  conclusions 
Parkin  and  Kermeen  mentioned  the  definite  need  for  further  work  along 
their  lines  in  order  to  explain  theoretically  the  experimentally  observed 
scale  effects.  Almost  no  experimental  work  along  the  lines  of  Parkin 
and  Kermeen  has  been  carried  out  since,  but  many  theoretical  models 
have  been  proposed  utilizing  the  interpretations  of  Parkin  and  Kermeen*  s 
experimental  observations. 

Review  of  Existing  Theoretical  Models:  In  general  theoretical  predic¬ 
tions  of  cavitation  inception  are  based  on  calculations  of  the  growth  of 
a  spherical  bubble  from  an  initial  radius  Rq  to  a  prescribed  final  radius. 
The  growth  of  the  bubble  is  assumed  to  be  determined  by  the  solution 
of  the  Rayleigh-Plesset  equation^ ^  to  be  considered  later.  The  pres¬ 
sure  distribution  around  the  headform  is  assumed  to  be  that  existing 
under  non-cavitating  conditions.  The  behavior  of  a  single  spherical 
bubble  is  assumed  to  characterize  the  problem. 

Parkin* 8  Theory;  One  of  the  first  attempts  to  model  the 

(24) 

process  of  cavitation  inception  theoretically  was  made  by  Parkin'  1 
in  1952.  In  his  formulation  Parkin  did  not  consider  the  effects  of 
boundary  layer  presence  on  the  model  surface.  It  was  supposed  that 
cavitation  is  initiated  from  small  nuclei  which  contain  air  or  water 
vapor,  or  both,  stabilized  on  small  solid  particles  in  the  liquid.  Based 
on  stability  arguments  and  experimental  observations  these  nuclei  are 

_4 

taken  to  be  equivalent  to  a  spherical  bubble  of  radius  Rq  =  2  X  10 
millimeters.  This  bubble  is  assumed  to  travel  along  the  model  surface 
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at  the  local  liquid  velocity;  during  this  travel  the  growth  of  the  bubble 
was  estimated  from  the  Rayleigh-Ples  cet  equation.  This  equation 
ignoring  the  viscous  and  temperature  effects  in  the  bubble  growth  is^^ 
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where  R  is  the  instantaneous  bubble  radius  and  a  dot  represents  the 
time  derivitive.  is  the  cavitation  number  based  on  vapor  pressure 
and  the  pressure  coefficient  C^(t)  describes  the  time  varying  liquid 
pressure  seen  by  the  bubble  as  it  travels  along  the  model  surface. 
Thus,  the  pressure  distribution  on  the  model  under  non-cavitating 
conditions  is  considered  to  be  known.  The  value  of  P^  is  the  gas 
pressure  inside  the  bubble  which  is  assumed  to  follow  ideal  gas  law 
giving  Pq=Pq  (Rq/R)  where,  Pq  and  Rq  are  initial  conditions  to  be 
prescribed. 

Parkin  assumed  the  surface  tension  term  to  be  equal  to  zero 
initially  and  to  reach  its  full  value  of  2T/R  when  R=  5Rq.  Thus,  his 
assumption  is  not  quite  equivalent  to  a  free  stream  nuclei  in  the  form 
of  a  gas  bubble  of  radius  Rq.  Inception  of  cavitation  was  assumed  to 
take  place  when  the  bubble  radius  reaches  a  value  of  1  millimeter  or 
when  the  final  radius  reaches  500  times  the  initial  radius.  Calculations 
carried  out  by  this  method  for  the  incipient  cavitation  number  on 
hemispherical  nosed  models  showed  poor  agreement  with  the  observed 
values.  Parkin  attributed  the  lack  of  agreement  between  theory  and 
experiment  to  the  neglect  of  the  effect  of  the  boundary  layer.  It  should 
be  mentioned  that  experiments  of  Parkin  and  Kermeen^^  showing  the 
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importance  of  boundary  layer  in  the  process  of  cavitation  inception 
were  done  after  Parkin  formulated  his  theory. 

(251 

Holl  and  Kornhauser's  Theory:  Holl  and  Komhauser'  ' 
employed  similar  ideas  as  Parkin's,  but  they  did  extensive  studies 
regarding  the  effect  of  initial  conditions,  stream  versus  surface  nuclei, 
and  thermodynamic  effects.  Thermodynamic  effects  were  considered 
in  determining  the  cavitation  inception  number  for  a  large  range  bulk 
fluid  temperatures.  This  required  the  simultaneous  solution  of  the 
Rayleigh-Plesset  equation  for  the  bubble  growth  and  the  energy  equation 
for  the  bubble  wall  temperature.  These  computations  were  made  as 
follows:  for  a  prescribed  Rq,  from  static  stability  analysis  there  exists  a 
liquid  pressure,  P^  ,  for  static  equilibrium.  The  corresponding  posi¬ 
tion  on  the  model  surface  at  which  the  local  pressure  is  equal  to  P^ 
is  taken  to  be  the  starting  point  of  bubble  history.  If  the  nucleus  was 
a  free  stream  nucleus,  it  was  assumed  to  travel  at  the  local  velocity  of 
the  liquid,  Vgj  in  particular,  if  boundary  layer  effectr  are  included  the 
velocity  of  the  bubble,  Vg,  was  taken  to  be  equal  to  Vg  within  the 
boundary  layer  at  a  distance  of  one  bubble  radius  from  the  model  sur¬ 
face.  If  the  nucleus  was  of  the  surface  type,  Vg  initially  was  taken  to 
be  zero.  Subsequent  values  of  Vg  were  determined  by  carrying  out  a 
trajectory  analysis  of  the  bubble  within  the  boundary  layer  accounting 
for  inertial  effects,  viscous  drag,  and  buoyancy  forces  arising  from  the 
pressure  distribution  on  the  body.  The  velocity  of  the  liquid,  Vg,  was 
again  taken  at  a  distance  of  one  bubble  radius  from  one  model  surface. 
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Cavitation  inception  was  said  to  occur  when  the  bubble  reaches 
a  maximum  diameter  of  0.01  in.  These  calculations  carried  out  for  a 
hemispherical  body  showed  that,  if  one  includes  thermodynamic  effects 
in  the  process  of  cavitation  inception,  the  value  of  a.  should  decrease 
considerably  with  increase  in  fluid  temperature.  Experiments,  on  the 
contrary,  showed  an  increase  in  with  increase  in  fluid  temperature. 
This  lack  of  agreement  was  attributed  to  possible  effects  of  convective 
heat  transfer  in  the  solution  of  energy  equation  neglected  in  their 
theoretical  formulation.  It  should  be  mentioned  that  without  any  physical 
basis  whatsoever  they  assumed  turbulent  transition  of  the  boundary  layer 
to  take  place  at  the  minimum  pressure  point. 

The  Van  der  Walle  Theory:  The  Van  der  Walle^2^  theory  is 
again  basically  the  same  as  that  of  Parkin's  except  that  he  considers 
in  addition  the  possibility  of  gaseous  growth,  vaporous  growth  of  both 
free  stream  and  surface  nuclei,  and  a  hypothesized  process  of  a  surface 
nuclei  being  generated  from  a  crevice.  He  concluded  that  the  surface 
nuclei  generation,  their  detachment  and  subsequent  behavior  will  be 
influenced  by  a  balance  of  local  pressure  gradient  forces,  inertial  drag 
forces  and  surface  adhesion  forces.  In  the  calculation  of  the  drag  force 
on  the  attached  nucleus,  the  local  flow  was  considered  to  be  laminar  or 
turbulent  and  the  local  velocity  within  the  boundary  layer  was  used  in  a 
similar  manner  to  that  of  Holl' s  theory  explained  earlier  but,  in  fact, 
his  work  was  first.  Thus  for  different  bodies  with  different  pressure 
gradients.  Van  der  Walle  obtained  various  inception  relations  and, 
accordingly,  scaling  factors.  Finally,  comparison  of  his  theoretical 
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calculations  of  on  a  hemispherically  nosed  body  showed  only  a 
qualitative  agreement  with  the  experimental  results.  Boundary  layer 
calculations  on  this  body  were  estimated  with  elementary  flat  plate 
boundary  layer  growth  relationships. 

The  Oshima  Theory:  Oshima's  theory^^  is  somewhat  different 
in  the  sense  that  he  attempts  to  formulate  strict  scaling  relationships 
based  on  similarity  ideas  and  for  this  reason  the  value  of  Oj  must  be 
known  at  a  reference  condition  from  experiments.  Oshima  assumed 
that  cavitation  occurs  at  geometrical  similar  positions  on  geometrical 
similar  bodies.  Dynamically  similar  growth  of  vapor  bubbles  as 
specified  by  the  Raylaigh-Plesset  equation  is  assumed  to  occur  between 
two  geometrically  similar  points.  Vapor  bubbles  are  assumed  to 
travel  at  a  distance  from  the  body  surface  which  is  a  constant  fraction  of 
the  local  boundary  layer  height.  Finally,  cavitation  is  assumed  to  occur 
when  the  bubble  diameter  is  equal  to  that  of  the  local  boundary  layer 
displacement  thickness,  an  assumption  based  on  Parkin  and  Ker- 
meen'sv  '  observations.  This  last  assumption  introduces  a  Reynold 
number  effect.  Boundary  layer  transition  was  assumed  to  occur  at  the 
point  of  instability  of  the  laminar  boundary  layer.  From  this,  the  flow 
at  the  position  of  inception  could  either  be  laminar  or  turbulent.  Oshima's 
results  compared  well  with  the  experiments  even  though  as  we  will  see 
from  the  present  work  that  his  last  assumption  is  based  on  wrong  inter¬ 
pretations  of  Parkin's  experiments. 


Johnson  and  Hsieh's  Theory:  The  previously  considered 
theoretical  models  place  main  importance  on  the  bubble  growth  which 
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takes  place  inside  the  boundary  layer.  But  considerations  of  the  growth 

of  free  stream  gas  bubbles  into  macroscopic  cavitation  as  observed 

by  Knapp^  can,  in  fact,  disregard  the  presence  of  the  boundary  layer 

unless  these  bubbles  somehow  get  entrained  into  it.  It  was  mentioned 

earlier  that  if  there  is  a  sufficient  number  of  free  stream  nuclei  in  the 

-2 

form  of  gas  bubbles  with  radii  larger  than  say  10  in.  present  in 

the  stream,  the  value  of  Cj  would  be  expected  to  be  equal  to  -  Cpmjn 

from  bubble  stability  considerations.  Johnson  and  Hsieh^®^  consider 

the  trajectories  of  these  free  stream  bubbles  with  the  flow  past  a  body. 

They  find  that  under  the  effects  of  the  pressure  gradient  there  is  a 

certain  "screening"  process  of  the  various  sized  bubbles  which  tends  to 

deflect  the  larger  bubbles  away  from  the  minimum  pressure  area  near 

the  body.  This  screening  process  is  speed  and  size  dependent  and  leads 

to  scaling  relations  again  based  on  the  idea  of  a  "critical"  radius,  but 

with  certain  assumed  free  stream  nuclei  distributions.  The  consistent 

trend  of  these  relations  is  for  an  increase  of  cr  with  increase  in  speed 

(21) 

and  body  size.  Scheibe  7  has  used  these  above  mentioned  ideas  to 
compute  Gj  with  more  realistic  free  stream  nuclei  distribution  and  his 
computations  agreed  only  qualitatively  with  experimentally  measured 
values. 


Scope  of  Present  Work:  From  previous  comments  it  is  clear  that,  both 
experimental  and  theoretical  investigations  have  neglected  any  detailed 
consideration  of  possible  real  fluid  effects  on  cavitation  inception. 
Especially  noteworthy  oversights  are  the  possibility  of  separation  and 
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of  any  consideration  for  the  unsteady  effects  of  transition  itself. 

Transient  pressures  are  known  to  influence  the  value  of  ^  in  separated 

flows  behind  disks^^,  rearward  facing  steps^^,  fully  developed  tur- 

(31)  (32) 

bulent  channel  flow'  'and  in  the  mixing  layer  of  a  submerged  jet 

(33) 

as  discussed  by  Arndt  and  Ippen  in  their  review 

Through  private  communication  J.  W.  Holl  of  the  Pennsylvania 

State  University  has  brought  to  our  attention  some  recent  work  by 
(34) 

A.  B.  Bailey'  'and  M.  V.  Casey.  In  their  work  they  have  shown  a 
good  correlation  of  measured  a ^  with  negative  value  of  the  pressure 
coefficient  at  the  predicted  position  of  laminar  separation  on  hydro¬ 
foils  at  moderate  angles  of  attack.  The  position  of  separation  was 
estimated  by  Thwaites  approximate  method.  But,  the  agreement 
mentioned  above  does  get  worse  with  increase  in  Reynolds  number 
and  increase  in  angle  of  attack.  As  will  be  discussed  later  similar 
calculations  during  early  portion  of  present  investigation  prompted 
us  to  take  up  the  subsequent  experimental  study. 

It  may  be  mentioned  here  that  there  are  several  possible  viscous 
flow  regimes  on  a  streamlined  body;  these  include,  (a)  laminar  boundary 
layer  flow  without  separation  or  transition,  (b)  laminar  boundary  layer 
flow  with  transition  but  without  any  separation,  (c)  same  as  (b)  but  with 
possible  turbulent  separation  after  transition,  and  (d)  flow  with  laminar 
boundary  layer  separation  with  the  possible  subsequent  transition  in  the 
separated  free  shear  layer.  It  is  expected  that  inception  would  be 
dependent  on  the  flow  regime .  F urther,  it  might  be  misleading  to  draw 
comparisons  of  values  of  o.  obtained  in  different  flow  regimes  without 
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this  knowledge.  The  present  work  is  concerned  with  identifying  the 
type  of  viscous  flow  regime  on  some  of  the  bodies  used  in  cavitation 
research  by  a  new  technique  and  further  to  determine  how  cavitation 
inception  may  be  influenced  by  these  real  fluid  flows. 

Our  interest  in  this  approach  was  prompted  by  some  preliminary 
boundary  layer  growth  calculations  made  for  the  I-  T.  T.  C.  test  body. 

It  was  found  that  the  position  of  inception  roughly  coincided  with 
the  predicted  position  of  the  laminar  separation.  It  was  also  found  that 
the  value  of  the  inception  cavitation  number  was  very  close  to  the  pres¬ 
sure  coefficient  at  this  point.  These  findings  stimulated  us  to  take  a 
new  look  at  the  interaction  of  viscous  effects  with  mechanism  of  cavi¬ 
tation  inception  on  streamlined  bodies  along  the  same  lines  as  Parkin's^^ 
original  work.  It  was  believed  that  earlier  methods  of  flow  visualization 
used  in  water  tunnel  research  would  be  inadequate  and  that  a  different 
method  of  flor  .sualization  would  be  needed.  The  schlieren  technique 
widely  used  in  the  field  of  compressible  aerodynamics  was  subsequently 
(and  successfully)  developed  for  this  purpose^\ 

Electrolysis  bubbles  as  a  means  of  flow  visualization  is  of  com¬ 
mon  use  in  slow  liquid  flows.  For  the  present  studies  it  was  thought 
that  this  principle  could  be  utilized  as  a  source  of  artificial  surface 
nuclei.  A  modest  investigation  of  the  use  of  this  technique  for  research 
in  cavitation  inception  seems  to  be  in  order. 
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2.  EXPERIMENTAL  METHODS 

2,1  Flow  Visualization 

Need  for  New  Method;  Normally  conventional  methods  of  flow 
visualization  such  as  dye  injection  and  oil  film  techniques  are  difficult 
to  use  in  water,  especially  at  high  velocities.  The  author  had  come  to 
realize  these  difficulties  in  his  earlier  work  with  flow  visualization  on 
hydrofoils^)  by  oil  film  technique.  During  some  preliminary  work  of 
the  present  studies,  the  method  of  dye  injection  was  used  to  show  the 
existence  of  a  laminar  separated  flow  on  a  two  inch  hemispherical  nose 
at  a  tunnel  velocity  of  about  5  feet  per  second,  but  at  a  velocity  of  about 
15  feet  per  second  the  observations  were  difficult  to  interpret  and  were 
confusing.  Due  to  these  difficulties  a  better  method  of  flow  visualization 
had  to  be  developed.  Many  methods  were  considered  including  inter¬ 
ferometry,  schlieren  technique  and  skin  friction  gages.  The  latter  were 
thought  to  be  feasible  but  expensive  to  obtain  spatial  resolution.  Inter¬ 
ferometry  and  the  schlieren  technique  both  require  a  density  gradient 
in  the  vicinity  of  the  viscous  boundary  layer.  It  was  though*  the  necessary 
gradients  could  be  created  by  heating  or  cooling  of  the  model  and  that 
the  resulting  thermal  boundary  layer  could  then  be  visualized  by  inter¬ 
ferometry  or  the  schlieren  technique.  Of  the  two,  the  latter  is  far 
simpler  and  has  a  long  history  of  use  in  compressible  aerodynamics. 

The  gross  properties  of  thermal  and  hydrodynamic  boundary  layers 
especially  when  visualizing  such  details  as  separation  and  transition 
are  very  similar.  Finer  details  such  as  relationship  between  the  growth 
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of  the  two  boundary  layers  are  dependent  on  the  Prandtl  number  and  the 
starting  position  of  heating  or  cooling. 

Preliminary  Work:  The  feasibility  of  the  schlieren  method  was  deter¬ 
mined  by  carrying  out  preliminary  experiments  in  the  Free  Surface 
Water  Tunnel  (FSWT).  The  working  section  of  this  tunnel  is  20  in.  wide 
and  8  ft.  long  with  a  water  depth  of  21  in.  ;  the  side  walls  of  the  tunnel 
are  made  of  transparent  plexiglass.  Further  details  of  the  tunnel  may 
be  found  in  Ref.  (36).  The  physical  principles  in  the  schlieren  method  of 
flow  visualization  are  well  known  and  may  be  found  in  Ref.  (37).  A 
schematic  drawing  of  the  optical  set-up  used  for  schlieren  purposes  in 
the  FSWT  is  shown  in  Fig.  3.  The  light  source  used  was  an  air  cooled 
steady  mercury  vapor  lamp  which  in  conjunction  with  a  gathering  lens 
and  slit  arrangement  as  shown  in  Fig.  3  formed  a  point  source.  The 
position  of  the  focusing  lens  and  the  film  plate  was  suet,  that  a  5:1 
magnification  could  be  obtained  with  adequate  illumination. 

The  first  experiments  were  on  a  two  inch  hemispherical  nose 
model  and  the  density  gradient  required  for  the  schlieren  effect  was 
accomplished  by  internal  circulation  of  hot  water  with  an  inlet  tempera¬ 
ture  of  130°F.  The  general  arrangement  and  the  position  of  the 
measurement  of  various  temperatures  to  obtain  heat  transfer  data  are 
shown  in  Fig.  4. 

Observations  made  with  the  continuous  light  source  showed  the 
existence  of  a  laminar  boundary  layer  separation  on  the  hemispherical 
nose  even  at  the  highest  tunnel  velocity  attainable  in  the  FSWT  (25  feet 
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per  second).  Photographs  obtained  with  the  help  of  a  mechanical  shut¬ 
ter  which  provided  exposure  time  of  1/100  sec.  did  not  show  any  of  the 
details  of  the  separated  free  shear  layer  or  its  subsequent  reattachment 
process,  if  the  latter  did  exist.  To  obtain  such  details  the  mercury 
vapor  lamp  was  replaced  by  a  spark  gap  source  which  had  an  exposure 
time  of  about  5  microseconds.  The  spark  gap  was  made  up  of  two  brass 
disks  each  of  0.5  in.  in  diameter  and  0.  14  in.  in  thickness  having  a 
minimum  clearance  of  about  0.2  in.  A  potential  of  10-15  kilovolts 
applied  to  the  disks  provided  a  bright  spark  of  the  time  duration  noted. 
The  spark  so  created  was  masked  to  provide  a  diffuse  source  light  of 
about  0.2  in.  X  0.14  in..  As  shown  in  Fig.  3  this  extended  source  was 
placed  at  the  focus  of  the  columnating  lens.  It  may  be  mentioned  that  a 
cut-off  plate  is  used  in  the  schlieren  apparatus  as  shown  in  Fig.  3.  This 
plate  had  to  be  adjusted  to  obtain  maximum  contrast  and  sensitivity  to 
provide  photographs  with  small  temperature  differences  needed  in  the 
present  experiment.  Shadowgraph  pictures  with  little  detail  could  be 
obtained  at  very  low  tunnel  speeds  by  removing  the  cut-off  plate  but, 
they  were  not  satisfactory  at  the  speeds  and  temperature  differences 
required  of  the  contemplated  experiment. 

The  photographs  obtained  with  this  schlieren  system  revealed 
clearly  the  laminar  boundary  layer  separation;  the  short  exposure 
time  resolved  some  of  the  features  of  the  subsequent  turbulent  re¬ 
attachment  process  at  these  tunnel  velocities  and  an  early  example  is 
provided  in  Fig.  11. 
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With  this  initial  success,  further  tests  were  carried  out  to 
estimate  heat  transfer  requirements  at  higher  tunnel  velocities  than 
used  in  the  FSWT.  With  the  help  of  an  iron-constantan  thermocouple 
imbedded  at  60  degrees  from  the  stagnation  point  on  the  model  surface 
as  shown  in  Fig.  4,  it  was  concluded  that  a  temperature  difference  of 
about  5°F  is  needed  to  produce  visible  schlieren  effects.  Heat  transfer 
requirements  to  produce  such  a  difference  were  estimated  by  measur¬ 
ing  the  temperature  drop  of  the  circulating  hot  water  as  a  function  of 
tunnel  velocity.  From  these  measurements  it  was  predicted  that  for  a 
two  in.  hemispherical  nose  with  an  afterbody  similar  to  that  shown  in 
Fig.  4  a  heat  input  of  about  3000  BTU  per  hour  would  be  required  at  a 
tunnel  velocity  of  60  feet  per  second.  This  in  fact  limited  the  velocity 
to  be  used  for  the  subsequent  experiments  in  the  High  Speed  Water 
Tunnel  (HSWT)  to  this  value.  This  heat  transfer  requirement  could  just 
be  met  without  great  difficulty  by  electrical  heating. 

To  test  qualitatively  the  effect  of  temperature  on  the  position  of 
boundary  layer  separation,  additional  schlieren  photographs  were  made 
with  circulation  of  cold  water  with  an  inlet  temperature  of  50  °F.  As 
will  be  discussed  in  a  later  section  there  are  no  significant  differences 
observed  for  the  conditions  of  the  subsequent  tests. 

It  should  be  mentioned  that  the  schlieren  technique  has  been 
previously  shown  to  work  in  water  by  Bland  and  Pelick^)  but  they  did 
not  subsequently  pursue  the  method. 


Final  Schlieren  Set-Up:  Success  of  the  preliminary  work  prompted  us 
to  try  the  method  in  the  HSWT.  To  provide  a  vibration  free  optical 
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system,  the  schlieren  set-up  was  mounted  on  a  platform  suspended 
from  the  ceiling  wall.  Optical  components  used  for  the  final  set-up 
were  the  same  as  that  used  for  the  preliminary  work,  but  with  some 
modifications  which  will  be  considered  next.  Basically,  the  schlieren 
apparatus  consisted  of  two  light  sources,  one  a  short  duration  spark 
gap  and  the  other,  a  steady  mercury  vapor  lamp.  The  steady  source 
was  utilized  for  visual  flow  observations  and  also  for  focusing  the 
optical  system.  A  schematic  diagram  of  the  apparatus  is  shown  in 
Fig.  5.  The  removable  mirror  shown  in  this  diagram  was  used  to  align 
the  spark  gap  system  with  the  aid  of  the  steady  light  source. 

It  may  be  noted  from  Fig.  5  that  the  windows  of  HSWT  are 
cylindrical  and  not  flat  as  in  the  FSWT.  This  introduces  severe  astig¬ 
matism  which  must  be  corrected  to  obtain  acceptable  schlieren  photo¬ 
graphs.  This  is  accomplished  by  placing  a  piano  convex  lens  on  each  of 
the  tunnel  side  windows  to  remove  first  order  astigmatism.  In  addition 
the  plexiglass  outer  surface  of  the  windows  was  not  optically  finished 
and  had  a  pronounced  "orange  peel"  surface  which  overwhelmed  the 
sought  for  schlieren  effect.  This  problem  was  overcome  to  a  large 
extent  by  introducing  a  thin  film  of  glycerine  between  the  correction 
lenses  and  the  window.  For  these  reasons  the  quality  of  the  photographs 
obtained  in  the  HSWT  were  not  as  good  as  that  obtained  in  the  FSWT, 
but  they  were  of  sufficient  quality  to  visualize  the  flow  at  all  velocities. 

A  photograph  of  the  apparatus  on  the  light  sources  side  is  given 
in  Fig.  6  and  Fig.  7  shows  the  camera  side. 
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The  required  heating  of  die  models  was  accomplished  by  instal¬ 
lation  of  a  Thunderbolt  kilowatt  electric  cartridge  heater  TB703 
(supplied  by  the  Vulcan  Electric  of  ITT)  within  the  model.  Details  are 
provided  in  Fig.  8. 

2.2  Test  Models 

Two  axisymmetric  models  widely  utilized  in  cavitation  inception 
studies  and  commonly  known  as  Hemispherical  nose  and  Swedish  head- 
form  (or  I.  T.  T.  C.  standard  headform)  were  used  in  the  present  investi¬ 
gations.  Details  of  the  test  bodies  are  shown  in  Figs.  9  and  10.  The 
Swedish  headform  presently  used  is  the  same  one  used  by  Acosta  and 
Hamaguchi  in  their  inception  studies  detailed  in  Ref.  (9).  The  test  bodies 
were  mounted  in  the  axisymmetric  test  section  of  the  HSWT^'  ^  by  a 
three  bladed  sting  support,  which  was  about  seven  diameters  downstream 
of  the  model  nose.  This  mounting  provided  a  rigid  support  for  the 
model  and  during  all  tests  no  vibration  or  movement  of  the  model  could 
be  discerned.  The  circular  cross  section  of  the  tunnel  was  14  inches 
in  diameter,  thus  the  models  occupied  less  than  two  percent  of  the 
through  flow  area.  The  approximate  angular  misalignment  between  the 
model  axis  and  tunnel  axis  was  measured  to  be  about  0.2 9  for  both  the 
test  bodies.  The  models  were  cleaned  with  taluene  before  and  after  each 
test  run. 


General 


>erimental  Procedures 


Phctographic  Methods:  Most  of  the  photographs  were  taken  under 
silhouette  conditions,  and  the  lens  and  film  plate  positions  were  such 
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that  a  magnification  of  5:1  was  obtained.  Single  exposure  photographs 
were  taken  with  the  spark  gap  which  had  an  exposure  time  of  about  5 
microseconds.  This  short  time  provided  freezing  action  of  the  physical 
processes  such  as  movement  of  cavitation  bubbles  and  growth  of  insta¬ 
bility  waves  on  the  separated  free  shear  layer  at  all  velocities,  thus 
producing  photographs  with  great  clarity.  Some  photographs  were  taken 
with  usual  arc  discharge  lamps  placed  on  the  camera  side  of  the  tunnel 
as  done  in  Ref.  (9).  High  speed  kodak  5  inches  X  7  inches  and  4  inches  x 
5  inches  Royal  x-pan  sheet  films  with  an  ASA  rating  of  1250  were  used 
for  single  exposure  photographs. 

Data  from  developed  negatives  concerning  positions  of  boundary 
layer  separation,  cavitation  separation  and  so  on  were  obtained  with  the 
help  of  a  negative  containing  nose  of  the  model  as  a  reference  point  and 
calibration  factor  obtained  by  photographing  a  scaled  negative  in  the 
place  of  the  model. 

Some  motion  pictures  of  the  cavitation  bubble  movements  were 
taken  with  a  Fastax  camera  at  4000  frames  per  second  and  few  at  5000 
frames  per  second.  Framing  rate  was  limited  by  the  available  intensity 
of  the  light  at  the  recording  film  position  to  produce  sufficient  contrast 
on  the  developed  negative.  The  light  source  used  was  a  steady  high 
pressure  Peck  110  mercury  vapor  lamp.  Motion  picture  studies  were 
recorded  with  16  mm  kodak  2475  Recording  Film  and  were  developed 
in  Accufine  solution  to  obtain  better  contrast.  Data  from  developed  film 
concerning  bubble  sizes  were  obtained  in  a  similar  manner  as  in  the 
case  of  single  exposure  negatives  mentioned  earlier. 
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Electrolysis  Bubble  Production:  Use  of  electrolysis  bubbles  for  flow 
visualization  has  been  widely  used  but,  in  present  studies  it  was  decided 
to  use  these  bubbles  as  source  of  artificial  surface  nuclei.  For  this 
purpose  the  whole  model  and  the  supporting  sting  were  used  as  two 
terminals.  A  non-conducting  adapter  was  utilized  to  electrically 
isolate  the  model  from  the  sting  assembly.  A  voltage  pulse  of  75  volts 
in  magnitude  and  approximately  1  second  in  duration  was  applied  between 
the  terminals  to  produce  the  electrolysis  bubbles.  This  produced  a 
total  current  of  15  amperes  in  the  circuit  and  a  current  density  of  about 
.2  amperes  per  square  inch  of  the  model  surface. 

Tunnel  Operation;  Pressure  measurements  for  the  computation  of  the 
free  stream  velocity  and  cavitation  number  were  made  with  mercury 
manometers.  As  there  is  a  boundary  layer  growth  along  the  walls  of 
the  working  section  of  the  tunnel,  it  is  necessary  to  calibrate  the  tunnel 
both  bare  and  with  a  model  installed  in  order  to  determine  the  static 
pressure  and  the  velocity  head  at  the  location  of  the  model.  Such  cali¬ 
brations  have  been  done  by  Acosta  and  Hamaguchi  and  their  results  are 
detailed  in  Appendix  A  of  the  Ref.  (9).  Working  formulae  for  the  compu¬ 
tation  of  tunnel  speed  and  cavitation  number  are  given  in  the  same 
Ref.  (9),  and  these  were  utilized  for  present  purposes  for  both  the 
models.  These  corrections  do  not  include  blockage  effects  and  may  be 
justified  from  the  computations  of  Hoyt^  who  has  shown  that  for  the 
Swedish  headform  when  the  model  occupies  about  two  percent  of  the 
through  flow  area  (as  is  the  case  for  present  experimental  conditions) 
the  blockage  effects  are  negligible. 
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Air  content  value  of  the  tunnel  water  was  held  fixed  at  a  value 
of  about  10  PPM  during  all  test  runs.  At  the  end  of  each  day  of  testing, 
the  air  content  value  was  measured  with  a  Van  Slyke  blood  gas  analyzer. 
At  this  particular  value  of  air  content  none  or  very  few  visible  macro¬ 
scopic  air  bubbles  were  seen  in  the  water  approaching  the  test  model 
under  all  conditions.  This  being  the  consequence  of  the  presence  of  a 
resorber  in  the  present  tunnel  facility.  Lower  air  content  would  have 
been  preferable  but  could  not  be  achieved  due  to  a  leak  in  the  main  seal 
of  the  pump. 

All  the  tests  were  conducted  at  room  temperature  and  the  cavi¬ 
tation  number  was  based  on  the  vapor  pressure  at  room  temperature. 

Method  of  Inception  Observations:  The  incipient  cavitation  conditions 
were  obtained  by  holding  the  tunnel  velocity  fixed  and  gradually  reducing 
the  tunnel  static  pressure  until  the  presence  of  steady  cavitation  near  or 
on  the  body  was  visually  observed  under  the  illumination  of  stroboscopic 
light.  At  this  point  all  the  manometer  readings  were  recorded  and  sin¬ 
gle  silhouette  photographs  were  taken  with  the  schlieren  apparatus.  By 
lowering  the  tunnel  static  pressure  further,  process  of  development  of 
cavitation  was  photographed  in  a  similar  manner. 

As  mentioned  earlier  none  or  very  few  macroscopic  bubbles 
were  observed  in  the  free  stream  approaching  the  test  model;  occa¬ 
sionally  these  cavitated  near  the  model,  but  this  type  of  cavitation  was 
avoided  and  not  recorded  during  present  studies.  In  the  terminology  of 
Parkin  and  Kermeen  on  Page  29  of  Ref.  (22),  present  incipient  conditions 
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may  be  termed  "boundary  layer  incipient  cavitation", 
were  made  to  record  desinent  conditions. 


No  attempts 
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3.  EXPERIMENTAL  RESULTS 

Flow  Visualization  Under  Fully  Wetted  Conditions;  Figure  11  shows  a 
schlieren  photograph  of  laminar  boundary-  layer  separation  bubble  at  a 

5 

Reynolds  number  of  10  on  a  tuo  inch  hemispherical  nose  mounted  in 
the  FSWT.  The  top  photograph  in  Fig.  12  shows  the  boundary  layer 
separation  when  the  surface  temperature  of  the  model  is  about  10 °F 
below  the  free  stream  water  temperature,  whereas  the  bottom  photo¬ 
graph  illustrates  the  same  phenomenon  when  the  surface  temperature  is 
about  10CF  above.  Both  photographs  are  at  a  Reynolds  number  of 
2.5  X  105  and  were  taken  in  the  same  facility.  As  mentioned  earlier 
the  difference  in  temperature  under  heated  conditions  was  measured 
directly,  but  tinder  cooled  conditions  it  is  estimated  from  heat  transfer 
measurements.  Schlieren  photographs  in  Fig.  13  were  taken  in  the 
HSWT  and  are  meant  to  show  the  behavior  of  laminar  separation  bubble 
with  changes  in  Reynolds  number  and  body  shape.  Also  shown  in  Fig.  13 
are  two  photographs  taken  without  heating  the  test  models  at  a  tunnel 
veloci'y  of  40  feet  per  second.  Variation  of  measured  laminar  boundary 
layer  separation  position  on  the  hemispherical  nose  with  tunnel  velocity 
is  shown  in  Fig.  14.  Similar  results  for  the  Swedish  headform  are  shown 
in  the  next  Fig.  15.  These  measurements  were  obtained  from  schlieren 
photographs  and  an  example  of  the  position  of  separation  used  is  marked 
on  the  photograph  (c)  of  Fig.  13  by  an  arrow  and  "S".  Also  included  in 
Fig,  14  are  results  obtained  from  tests  conducted  in  the  FSWT  on  a 
different  hemispherical  nose  than  used  in  the  HSWT.  A  plot  of  non- 
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dimensional  axial  distance  from  the  point  of  separation  to  the  dis¬ 
appearance  of  separated  white  line  in  the  photographs  of  Fig.  13  as 
function  of  tunnel  velocity  is  shown  in  Fig.  16.  An  example  of  point  of 
disappearance  used  is  marked  on  the  photograph  (i)  of  Fig.  13  by  an 
arrow  and  "T".  A  measure  of  vertical  extent  of  the  laminar  separation 
bubble  was  obtained  from  the  measured  value  of  the  maximum  height 
from  the  body  surface  to  the  separated  white  line  in  the  photographs  of 
Fig.  13.  An  example  of  this  is  provided  in  the  photograph  (I)  of  Fig.  13 
and  is  marked  by  "H^".  Variation  of  non-dimensional  Hm  with  tunnel 
velocity  is  shown  in  Fig.  17  for  the  hemispherical  nose  and  also  for  the 
Swedish  headform. 

Cavitation  Inception  Studies:  Figure  18  provides  a  sequence  of 
photographs  illustrating  the  development  of  cavitation  in  the  separated 
layer  of  the  hemispherical  nose  at  a  velocity  of  40  feet  per  second. 
Partial  development  at  60  feet  per  second  is  shown  in  the  next  Fig.  19. 
The  dark  patches  above  the  model  surface  in  the  photographs  of  Figs. 

18  and  19  represent  the  cavitating  areas.  Figure  20  shows  two  different 
kinds  of  incipient  cavitation  data  for  the  hemispherical  nose.  The  un¬ 
filled  circles  in  Fig.  20  correspond  to  the  incipient  conditions  as  shown 
by  the  two  top  photographs  of  Figs.  18  and  19.  The  filled  triangles 
correspond  to  the  incipient  conditions  as  shown  by  the  third  photograph 
of  Fig.  18.  Comparison  of  present  data  with  that  obtained  by  Kermeetf^ 
is  shown  in  Fig.  20.  Variation  of  incipient  cavitation  number  for  the 
Swedish  headform  with  tunnel  velocity  is  shown  in  Fig.  21.  Comparison 
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of  average  data  of  Fig.  21  with  that  obtained  in  Ref.  (9)  is  shown  in  the 
next  Fig.  22. 

Electrolysis  Studies:  Figure  23  shows  a  series  of  photographs 
illustrating  the  effect  of  producing  artificial  nuclei  by  electrolysis  on 
the  hemispherical  nose  at  the  same  velocity  of  28  feet  per  second,  but 
at  various  cavitation  numbers.  The  photograph  marked  "natural 
cavitation"  is  taken  without  the  production  of  electrolysis  bubbles.  A 
similar  set  of  photographs  for  the  Swedish  headform  is  shown  in  Fig.  24 
These  photographs  were  taken  with  the  help  of  two  arc  dishcarge  lamps 
located  on  the  same  side  of  the  tunnel  as  camera  and  a  silhouette  photo¬ 
graph  was  superimposed  on  the  same  negative.  The  latter  was  done 
with  the  spark  gap  under  fully  wetted  flow  conditions.  Next  series  of 
photographs  in  Fig.  25  were  taken  at  about  the  same  cavitation  number, 
but,  at  different  velocities. 

Motion  Picture  Studies:  Some  interesting  cavitation  bubble  motions  in 
the  separated  layer  of  the  hemispherical  nose  are  shown  in  Figs.  26(a) 
and  26(b).  Both  series  of  photographs  were  taken  at  a  velocity  of  40 
feet  per  second  and  cavitation  number  of  0.74.  At  the  top  of  the  sequence 
is  a  schlieren  photograph  showing  the  boundary  layer  separation  at  40 
feet  per  second.  These  photographs  were  obtained  with  the  spark  source 
and  by  turning  the  film  manually  frame  by  frame.  Sequence  of  dis¬ 
appearance  of  cavitation  on  the  Swedish  headform  at  40  feet  per  second 
and  approximate  cavitation  number  of  0.46  is  shown  in  Fig.  27.  The 
schlieren  photograph  at  the  top  of  the  sequence  is  at  30  feet  per  second. 


4.  DISCUSSION  OF  FULLY  WETTED  FLOW  OBSERVATIONS 


Position  of  Boundary  Layer  Separation:  Effectiveness  of  the  schlieren 
method  as  flow  visualization  technique  is  clearly  illustrated  in  the 
photograph  of  Fig.  11.  It  shows  the  position  of  boundary  layer  separa¬ 
tion,  and  also  the  process  of  reattachment  of  the  separated  free  shear 
layer  to  form  a  laminar  separated  bubble. 

Measurements  shown  in  Fig.  14  illustrate  that  the  position  of 
boundary  layer  separation  did  not  alter  with  change  in  surface  tempera¬ 
ture  of  the  model  for  all  test  velocities  in  the  FSWT.  Similar  results 
were  observed  concerning  the  extent  of  separated  bubble  as  may  be 
seen  from  the  photographs  of  Fig.  12.  The  observed  difference  in  the 
break  up  mechanism  of  the  free  shear  layer  is  not  due  to  the  temperature 
difference,  but  is  an  inherently  unsteady  property  of  the  separated  layer. 
Thus  within  the  range  of  temperatures  used  to  create  the  schlieren 
effects,  the  position  or  the  extent  of  laminar  separated  bubble  was  not 
influenced  by  heating  or  cooling  of  the  model. 

As  may  be  seen  from  the  set  of  photographs  in  Fig.  13,  both 

hemispherical  nose  and  Swedish  headform  possessed  laminar  separation 

bubble  at  all  the  test  velocities  in  the  HSWT.  The  highest  test  velocity 

5 

corresponds  to  a  Re^  of  9.06  X 10  for  the  hemispherical  nose  and  ReD  of 
8.78  X  105  for  the  Swedish  headform.  To  help  distinguish  the  separated 
white  line  in  the  photographs  of  Fig.  13,  two  photographs  taken  without 
heating  of  the  model  are  included  at  the  top  of  the  sequences. 

From  Fig.  14,  it  may  be  seen  that  the  position  of  laminar 
boundary  layer  separation  on  the  hemispherical  nose  mounted  in  both  the 
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HSWT  and  FSWT  is  about  87*  from  the  stagnation  point.  To  avoid  free 
surface  effects  in  the  FSWT,  the  model  was  submerged  about  6  diameters 
below  the  surface.  It  should  be  noted  that  two  different  models  were 
utilized  in  the  two  different  tunnels  whose  characteristics  (such  as  free 
stream  turbulence  levels)  are  not  completely  known  and  may  be  different. 
In  both  tunnel  facilities  the  position  of  separation  on  the  hemispherical 
nose  did  not  change  appreciably  with  increase  in  velocity,  the  change 
being  within  the  error  of  experimental  measurements.  As  seen  from 
Fig.  15,  for  Swedish  headform  s/D  at  separation  is  equal  to  0.82  at  a 
velocity  of  20  feet  per  second  and  increases  to  0.84  at  a  velocity  of 
60  feet  per  second. 

Length  of  Separated  Bubble:  In  almost  all  the  photographs  of  Fig.  13, 
the  separated  white  line  begins  to  form  into  a  wavy  line  and  subsequently 
either  disappears  or  becomes  very  faint.  Since  the  schlieren  effects 
are  a  measure  of  the  density  gradients  (in  this  case  induced  by  tempera¬ 
ture  gradients)  disappearance  of  the  separated  white  line  may  be 
interpreted  as  vanishing  of  these  gradients  due  to  turbulent  mixing.  The 
distance  from  the  position  of  separation  to  the  point  of  disappearance  of 
the  separated  white  line  may  be  used  as  a  convenient  measure  of  the 
length  of  the  separated  bubble.  Variation  of  such  distance  with  velocity 
for  both  the  test  models  is  shown  in  Fig.  16.  As  expected  the  bubble 
length  decreases  with  increase  in  velocity  for  both  the  models.  Physi¬ 
cally  for  both  the  models  these  bubble  lengths  were  quite  small,  being 
only  12  percent  of  the  diameter  at  20  feet  per  second  (Rep£3.02  X  105) 
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and  about  4  percent  at  60  feet  per  second  (B.e^=  9  X 10  ). 

Some  Comparisons  With  Wind  Tunnel  Data:  Length  of  separated  bubble  - 

Lr  for  both  the  models  in  terms  of  appropriate  non-dimensional  variables- 

Rea  ,  L  /  9  is  shown  in  Fig.  28.  The  value  of  the  momentum  thickness 
08  r  s 

at  separation-9g  was  estimated  by  Thwaite’s  method  as  outlined  in 
Appendix  A.  Also  shown  in  Fig.  28  is  data  obtained  by  Gaster^0^  in  his 
wind  tunnel  studies  concerning  laminar  separation  bubbles  which  were 
created  on  a  flat  plate  with  the  help  of  an  auxilary  airfoil.  Distance 
from  the  position  of  separation  to  a  downstream  point  where  velocities 
close  to  the  wall  become  significantly  larger  (say  a  tenth  of  the  free  stream 
velocity)  was  used  as  a  measure  of  length  of  the  separated  bubble.  Gas- 
ter  refers  to  this  distance  as  "dead  air  region"  and  measurements  of 
were  obtained  with  the  help  of  detailed  velocity  and  pressure  measure¬ 
ments.  Favorable  comparison  with  Gaster's  data  indicates  that  present 
basis  and  method  of  obtaining  the  length  of  separated  bubble  is  quite 
adequate.  It  can  also  be  inferred  from  Fig.  28  that  Lr/  9g  is  slowly 

varying,  but  primarily  a  function  of  Re  0  . 

°s 

Height  of  the  Separated  Bubble;  As  seen  from  Fig,  13,  on  present 
bodies  the  height  of  the  separated  bubble  from  the  model  surface 
essentially  starts  with  a  value  of  zero  and  gradually  increases  to  a 
maximum  value  -Hm  at  the  end  of  the  separated  region.  From  Fig.  17, 
it  is  seen  that  Hm  decreases  with  increase  in  velocity  for  both  the 
models.  The  largest  physical  value  of  Hm  indicated  by  the  schlieren 
photographs  is  only  about  one  hundredth  of  the  body  diameter. 


X&g&t Srw***-«-*« 
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Measurements  shown  in  Fig.  17  are  within  25  percent  accuracy.  A 
replot  of  Fig.  17  in  terms  of  appropriate  non-dimensional  variables 
-  Re  q  ,  H^/fig  is  shown  in  Fig.  29.  At  a  large  Reynolds  number  of 

8  c 

about  8.8  X  10  ,  on  the  Swedish  headform  the  maximum  separated 
height  is  only  twice  the  value  of  estimated  momentum  thickness  at 
separation.  The  pertinent  point  to  note  here  is  that  at  the  same 
Reynolds  number  the  observed  maximum  separated  height  on  the 
hemispherical  nose  is  two  to  three  times  larger  than  on  the  Swedish 
headform;  this  being  the  case  for  all  the  test  Reynolds  numbers. 

A  Remark  on  Measurement  of  Data  From  Schlieren  Photographs:  In 
present  experiments,  height  of  the  separated  white  line  in  the  photo¬ 
graphs  of  Fig.  13  is  not  a  true  measure  of  the  actual  separated  height. 
The  exact  relationship  between  the  two  is  a  function  of  the  optical  set¬ 
up.  Similar  problems  are  encountered  in  locating  the  actual  shock 
position  from  schlieren  photographs  in  the  field  of  compressible  aero- 
dynamics'  .  These  problems  arise  only  with  measurements  which 
are  normal  to  the  plane  of  the  knife  edge.  The  quantities  measured 
which  are  in  parallel  with  the  plane  of  the  knife  edge  reflect  a  true 
measurement.  Thus,  the  results  shown  in  Figs.  14  and  15  indicate  the 
true  position  of  the  boundary  layer  separation.  Similar  arguments  hold 
for  the  results  in  Fig.  16,  but  only  comparative  studies  regarding  the 
separated  height  shown  in  Fig.  17  may  be  made. 


Comparison  With  Approximate  Computations:  As  outlined  in  Appendix 
A,  some  theoretical  estimates  of  the  position  of  boundary  layer  separation 


-38- 


using  Thwaites  approximate  method  were  made.  For  the  hemispherical 
nose  the  predicted  value  was  found  to  be  87*  and  agrees  very  well  with 
the  observed  value  of  87*  in  both  HSWT  and  FSWT.  Predicted  value 


of  s/D  at  separation  for  the  Swedish  headform  was  found  to  be  0.82  for 

all  Reynolds  numbers  whereas,  the  observed  value  changes  from  0.82  to 

0.84  with  increase  in  Rejj  from  2.7  X  10®  to  8.78  X  10®. 

It  is  well  known  that  at  high  enough  Reynolds  number,  attached 

laminar  boundary  layers  possess  a  point  of  instability  after  which  the 

outside  disturbances  grow  exponentially  in  amplitude  until  transition  to 

turbulence  sets  in.  A  semi-empirical  method  developed  by  Smith  and 
(411 

his  associates'  7  enables  one  to  estimate  the  approximate  position  of 
transition  and  use  of  this  method  to  predict  transition  on  present  bodies 
is  detailed  in  Appendix  B.  The  Reynolds  number  at  which  transition  is 
to  occur  at  the  position  of  laminar  boundary  layer  separation  may  be 


referred  to  as  a  "critical  Reynolds  number"  designated  by  "Re^  ". 


Presumably  for  ReD 
is  possible.  Roshko 


>  Ren  ,  no  laminar  boundary  layer  separation 
^crit 

(421 

'  7  refers  to  this  Reynolds  number  as  "second 


critical  Reynolds  number"  and  finds  it  to  be  about  5X10 


6 


for  the  case 


of  a  cylinder  in  rectilinear  flow.  Present  computations  from  Appendix 
B  indicate  that  R«DCrit^  5  X  10®  for  the  hemispherical  nose  and  that 
ReDcrit"**®  x  ^  *or  headform.  This  is  in  accord  with 


present  experimental  observations  that  the  boundary  layer  is  still 

,  5 

laminar  at  separation  at  a  ReD  =  9  X10  for  the  hemispherical  nose  and  at 
•  5 

a  Rej)=  8.7x10  for  the  Swedish  headform.  During  the  present  experi¬ 


ments  it  was  not  possible  to  attain  critical  Reynolds  number  for  either 
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of  the  models,  however  it  is  noteworthy  to  point  out  that  the  length  of 
the  separated  bubble  decreased  with  increase  in  Re^  and  length  was 
only  4  percent  of  the  model  diameter  at  Rej}  =  8.7  X  10^  for  the  Swedish 
headform. 

From  pressure  measurements  on  the  hemispherical  nose 
Rouse^^  indicates  that  &®Dcrit=  2  X  10^,  whereas  from  present 
experiments  R®DCrit  *8  *oun^  to  be  greater  than  9x10  .  It  is  believed 
that  this  discrepancy  has  resulted  from  the  use  of  pressure  taps  with 
a  spacing  of  12.5  percent  of  the  body  diameter  by  Rouse  in  the  measure¬ 
ment  of  the  pressure  distribution.  As  may  be  seen  from  Fig.  16,  the 
method  of  pressure  distribution  measurement  by  Rouse  would  not 

e 

indicate  the  possibility  of  laminar  separation  for  ReD  >  10  . 
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5.  DISCUSSION  OF  INCEPTION  STUDIES 

5.1  Hemispherical  Nose 

Visual  Observations:  Onset  of  cavitation  took  place  in  the  form  of  a 

ring  of  macroscopic  bubbles  around  the  nose  of  the  model.  This  type  of 

inception  may  be  termed  "steady  incipient  bubble  cavitation"  following  Ker- 
(4) 

meen's'  '  description  of  desinent  cavitation.  Close  observation 
indicated  the  presence  of  microscopic  bubbles  ahead  of  the  macro¬ 
scopic  ones  confirming  similar  observations  by  Parkin  and  Kermeen^’^ 
on  the  same  type  of  test  body.  As  tunnel  pressure  was  lowered  further, 
cavitation  bubbles  became  more  profuse  until  a  clear  macroscopic 
attached  cavity  was  formed  on  the  model  surface.  This  state  of 
cavitation  is  termed  "band  inception".  The  small  microbubbles 
characterizing  the  incipient  bubble  cavitation  disappear  although  at 

the  terminus  of  the  attached  cavity  such  bubbles  are  again  seen.  Ker- 
(4) 

meen'  '  refers  to  this  form  of  cavitation  as  "steady  attached  cavitation". 
The  development  just  described  was  found  to  take  place  for  all  the  test 

velocities  used  between  25  to  60  feet  per  second.  This  corresponds  to 

5  5 

Rejj  between  3.78  X 10  to  9.06x10  . 

Relationship  Between  Boundary  Layer  Separation  and  Inception:  With 
the  help  of  the  schlieren  method  it  was  possible  to  observe  the  inception 
process  and  the  nature  of  real  fluid  flow  past  the  model  simultaneously. 

As  can  be  seen  from  the  photographs  of  Figs.  18  and  19,  the  entire 
process  of  inception  takes  place  within  the  separated  region.  From 
photographs  18(a)  and  19(d),  it  is  evident  that  the  first  visible 


macroscopic  cavitation  occurs  in  the  transition  region  of  the  separated 
free  shear  layer  (as  mentioned  earlier  disappearance  of  the  separated 
white  line  is  taken  to  be  beginning  of  the  transition).  As  illustrated  in 
the  photographs  18(b)  and  19(e),  reduction  in  a  produces  increasing 
amounts  of  macroscopic  cavitation;  at  the  same  time  the  presence  of 
microscopic  bubbles  starting  from  the  position  of  laminar  separation 
becomes  apparent.  Further  reduction  in  c  results  in  the  formation  of 
a  single  macroscopic  attached  cavity  as  shown  in  the  photograph  18(c) 
which  also  shows  the  disappearance  of  the  incipient  bubble  cavitation 
observed  in  the  earlier  photographs.  It  may  be  noted  here  that  the 
position  of  laminar  separation  has  been  shifted  upstream  by  the  presence 
of  the  cavitation. 

It  should  be  noted  from  the  photographs  18(a)  and  19(d)  that  the 
diameters  of  the  macroscopic  bubbles  is  approximately  equal  to  the 
maximum  separated  height  -  H^.  Measurement  of  bubble  sizes  at  40 
feet  per  second  indicated  that  the  macroscopic  bubbles  were  of  the 
order  of  0.01  in.  in  diameter  and  the  microscopic  bubbles  were  of  the 
order  of  0.0025  in.  in  diameter.  These  values  are  in  agreement  with 
the  previous  measurements  of  Parkin  and  Kermeen^' 

From  the  photograph  in  Fig.  18(a),  it  is  important  to  note  that 
inception  of  cavitation  does  not  take  place  in  the  neighborhood  of  the 
minimum  pressure  point  which  is  expected  to  lie  upstream  of  the  posi¬ 
tion  of  laminar  separation. 

Further  Details  From  Motion  Picture  Studies:  To  get  further  detailed 
information  about  the  inception  process  some  motion  pictures  of 
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cavitation  in  the  separated  layer  were  taken  at  4000  frames  per  second, 

and  some  of  these  are  shown  in  the  sequential  photographs  of  Figs. 

26(a)  and  26(b).  These  photographs  were  taken  at  a  tunnel  velocity  of 

40  feet  per  second  and  a  a  of  0.74.  Also  included  at  the  top  of  the 

sequences  are  schlieren  photographs  taken  at  the  same  velocity  and  to 

the  same  scale.  From  frame-by-frame  study  of  these  motion  pictures 

it  could  be  deduced  that  both  macroscopic  and  microscopic  cavitation 

bubbles  orginated  in  the  transition  region  of  the  separated  layer;  the 

larger  macroscopic  bubbles  were  detached  from  the  model  surface  to 

be  collapsed  downstream  and  the  smaller  microscopic  bubbles  before 

they  could  grow  further  were  entrained  in  the  reverse  flow  of  the 

separated  region  and  traveled  upstream  towards  the  position  of  the 

separation  at  about  a  velocity  of  one  ft/sec.  as  compared  to  free  stream 

velocity  of  40  ft/ sec.  At  present  framing  rates  it  was  not  possible  to 

(22) 

track  the  bubbles  further.  Parkin'  '  however  observed  velocities  of 
these  small  bubbles  once  in  the  free  shear  layer  to  be  between  12  to  26 
ft/sec. 

These  events  were  the  ones  most  often  observed  but  there  were 
others  with  a  different  history.  A  few  bubbles  were  observed  to  grow  in 
the  separated  layer  at  a  fixed  position  on  the  body  surface  as  also  ob- 
served  by  Parkin  and  Kermeen'  if  the  growth  took  place  near  the 
position  of  separation  the  bubbles  became  entrained  into  the  free  shear 
layer  after  growing  to  a  size  equal  to  that  of  the  local  separated  height. 
If  the  growth  took  place  downstream  from  the  position  of  sep^ratfon  the 
bubble  motion  was  similar  to  that  described  initially  for  the  bubbles 
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which  started  from  the  transition  region.  Some  bubbles  were  seen  to 
circulate  twice  in  the  separated  region  before  they  were  entrained  in 
the  free  shear  layer  and  others  were  seen  to  start  from  the  transition 
region  to  be  transported  upstream  for  a  distance  after  which  they  were 
swept  downstream  although  still  within  the  separated  region.  An  example  of 
this  is  shown  in  the  photographs  of  Fig.  26(b).  The  important  point  to  note 
from  these  complex  bubble  motions  is  that  residence  time  of  the  bubbles 
in  the  separated  region  was  estimated  to  be  from  1  to  5  milliseconds 
and  during  this  period  they  did  not  show  any  appreciable  growth.  The 
bubbles  which  grew  at  a  fixed  position  on  the  model  surface  also  had  a 
growth  time  of  the  order  of  few  milliseconds  as  also  observed  by  Parkin 
and  Kermeen  but  the  growth  of  bubbles  in  the  transition  region  was  quite 
rapid  and  could  not  oe  photographed  in  any  detail  during  present  studies. 
Kermeen'  indicates  the  growth  time  of  these  latter  bubbles  to  be 
shorter  than  even  a  tenth  of  a  millisecond. 

Interpretation  of  Parkin's  Observations:  From  the  similarity  of  the 

(22) 

present  observations  and  those  formerly  carried  out  by  Parkin'  '  we 
may  make  the  following  important  points 

(a)  microscopic  bubble  growth  fixed  on  the  body  surface 
as  observed  in  Parkin’s  experiments  occurred  in  the 
laminar  separated  zone  and  not  submerged  in  the  attached 
boundary  layer  as  reported  by  Parkin, 

(b)  stabilization  of  these  bubbles  occurred  since  they 
were  in  the  separated  zone,  a  region  of  very  low  velocities 
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and  not  due  to  pressure  gradient  forces  as  thought  of 

(25) 

by  Parkin  and  subsequently  considered  by  Holl'  '  and 
Van  der  Walle^^  in  inception  scaling  theories, 

(c)  upstream  motion  of  these  bubbles  as  observed  by 
Parkin  were  due  to  the  reverse  motion  of  the  fluid  in  the 
separated  zone  and  not  due  to  adverse  pressure  gradient 
effects  as  thought, 

(22) 

(d)  it  is  mentioned  in  Parkin's  repdrt'  ’  that  the 
macroscopic  bubbles  are  sustained  by  the  microscopic 
bubbles  upstream,  but  present  experiments  indicate  that 

the  macroscopic  and  microscopic  bubbles  grew  independently 
of  each  other, 

(e)  after  a  clear  attached  cavity  was  achieved.  Parkin 
observed  the  suppression  of  the  microscopic  bubble 
growth.  It  is  believed  that  this  is  due  to  the  filling  of  the 
fully  wetted  separated  layer  with  the  cavity  leaving  only  a 
small  distance  of  about  .02  in.  between  the  position  of  viscous 
separation  and  the  beginning  of  the  cavity  which  may  be  seen 
in  the  photograph  (c)  of  Fig.  18.  Upstream  secondary  bubble 
motion  observed  by  Parkin  and  shown  in  the  Fig.  14  of  Ref. 

(22)  is,  in  fact,  within  this  small  region. 

It  should  be  mentioned  that  primarily  the  method  of  flow 
visualization  differed  between  the  present  experiments  and  that  con¬ 
ducted  earlier.  The  dye  streak  method  used  by  Parkin  was  not  adequate 
to  show  the  existence  of  a  small  laminar  separation  bubble;  instead  this 
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region  was  thought  to  be  the  transition  region  of  the  attached  boundary 
layer. 

Cavitation  Inception  Measurements:  From  Fig.  20,  increases  with 

increase  in  velocity  for  both  bubble  and  band  types  of  inception.  Since 

inception  was  seen  to  occur  downstream  of  the  position  of  laminar 

separation,  naturally  the  important  pressure  coefficient  to  compare  o^ 

with  would  be  the  negative  value  of  the  pressure  coefficient  at  separation, 

-  c  .  The  value  of  -  c  =  .63  shown  in  Fig.  20  is  obtained  with  the  help 
Ps  Ps 

(4) 

of  pressure  distribution  measurements  made  by  Kermeen'  '  and  the 
known  value  of  the  position  of  separation  in  the  HSWT  shown  in  Fig.  14. 
Also  shown  in  Fig.  20  is  the  value  of  -  Cpmjn=  .74. 

It  should  also  be  noted  that  0^  for  incipient  bubble  cavitation  at 
all  test  velocities  is  greater  than  -  Cp#,  indicating  that  cavitation  took 
place  at  local  static  pressures  higher  than  the  vapor  pressure.  This  is 
suggestive  that  either  gaseous  cavitation  was  observed  or  the  strong 
pressure  fluctuations  in  the  transition  region  were  responsible  for 
vaporous  cavitation.  From  motion  picture  studies  one  may  rule  out  the 
possibility  of  steady  gaseous  diffusion  as  the  mechanism  for  the  bubble 
growth,  since  the  residence  time  of  the  cavitation  bubbles  was  relatively 
short  and  this  point  has  been  discussed  by  Parkin'  '  in  some  detail. 

It  seems  possible  that  a  combination  of  convective  gaseous  diffusion 
rectified  gaseous  diffusion'  '  and  vaporous  cavitation  could  be  respon¬ 
sible  for  the  growth  of  macroscopic  bubbles  in  the  transition  region. 
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For  vaporous  cavitation  the  pressure  fluctuations  would  have  to 
be  of  the  order  of  12  percent  of  the  dynamic  head  at  30  feet  per  second 
and  25  percent  of  the  dynamic  head  at  60  feet  per  second.  Since  no 
pressure  fluctuation  measurements  under  the  present  type  of  flow 
situations  are  available,  further  discussion  on  the  origination  of 
cavitation  bubbles  in  the  transition  region  would  be  of  speculative 
nature. 

It  may  be  noted  from  Fig.  20  that  band  type  of  inception  also 

occurs  for  a.>  -c  indicating  that  this  type  of  inception  is  definitely 
1  Ps 

of  gaseous  typo. 

Comparison  with  Kermeen's  Data:  From  Fig.  20,  Kermeen's'  '  data 
of  desinent  cavitation  is  in  poor  agreement  with  present  values  of  a. 
for  incipient  bubble  cavitation,  but  is  in  good  agreement  with  present 
values  of  Oj  for  incipient  band  cavitation.  From  Kermeen's  report  it 
is  known  that  the  desinent  cavitation  was  of  the  bubble  type,  but  the 
state  of  cavitation  could  either  be  that  represented  by  Fig.  18(a)  or 
Fig.  18(b).  In  fact,  in  Fig.  20  the  uppermost  value  of  the  unfilled 
circle  represents  the  upper  bound  of  a.  for  incipient  bubble  cavitation, 
whereas  the  triangular  filled  values  represent  the  lower  bound.  A 
better  comparison  with  Kermeen's  data  will  be  made  later  in  Section  6 
and  is  shown  in  Fig.  30. 

Before  any  comparison  with  Kermeen's  results  is  made,  it 
should  be  made  certain  that  the  position  of  laminar  separation  did  not 
differ  appreciably  in  the  two  cases.  Of  course,  there  is  no  direct 
evidence  of  the  above  point  but  there  is  some  indirect  evidence.  In 
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Fig.  24  of  Ref.  (4)  a  photograph  is  shown  with  macroscopic  cavitation 

and  also  microscopic  cavitation  ahead  of  it;  the  starting  position  of  the 

latter  is  mentioned  to  be  88.7s  from  the  stagnation  point.  The  present 

studies  show  (see  Fig.  18(b))  the  same  type  of  cavitation  and  the  starting 

position  of  the  microscopic  bubbles  is  seen  roughly  to  coincide  with  the 

position  of  the  boundary  layer  separation.  Thus,  with  this  observation 

in  mind,  the  position  of  the  boundary  layer  separation  in  Kermeen's 

work  must  be  taken  to  be  about  88*  as  compared  to  87*  observed  on  the 

present  model.  In  this  event  -c  in  Kermeen's  case  would  then  be 

ps 

about  0.6. 


5.  2  Swedish  Headform 

Observations  of  Inception:  Inception  was  always  sudden  and  in  the  form 
of  a  fully  developed  cavity  or  was  of  band  type.  Typical  photographs  of 
inception  are  shown  in  (d)  of  Fig.  24  and  also  (d)  of  Fig.  25.  An  impor¬ 
tant  point  to  note  is  that  inception  takes  place  in  the  neighborhood  of  the 
position  of  the  boundary  layer  separation  (s/D  =  ,83)  which  lies  con¬ 
siderably  downstream  of  the  position  of  minimum  pressure  point^^^ 

(s/D  =  .62).  This  was  the  case  for  all  the  test  velocities  which  ranged 
from  25  to  60  feet  per  second,  corresponding  to  Re^  of  3.66  X 10  to 
8.78x10  .  Motion  pictures  taken  at  4000  frames  per  second  of  cavita¬ 
tion  in  the  separated  layer  produced  only  sketchy  details  of  the  inception 
process.  A  large  macroscopic  bubble  would  suddenly  appear  in  the 
separated  region  on  a  frame  of  the  motion  picture  film  and  by  the  next 
frame  a  fully  developed  cavity  would  be  formed.  However,  some 
details  of  the  desinent  cavitation  were  observed  as  shown  in  the  photo¬ 
graphs  of  Fig.  27  taken  at  a  velocity  of  40  ft/sec  and  o  =  .46.  It  may  be 
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seen  that  last  traces  of  cavitation  exist  in  the  transition  region  of  the 
separated  free  shear  layer.  It  would  be  too  premature  to  conclude 
from  this  that  inception  also  took  place  in  the  transition  region  of  the 
separated  free  shear  layer  even  though  on  the  hemispherical  nose 
inception  and  desinence  were  observed  to  take  place  in  the  same  region 
of  the  flow. 

Cavitation  Inception  Measurements:  From  Fig.  21,  except  at  a  velocity 
of  30  feet  per  second  the  inception  indexes  reproduce  quite  well.  At 
30  feet  per  second  considerable  scatter  exists.  Comparison  of  the 
average  data  of  Fig.  21  with  that  obtained  by  Acosta'  7  on  the  same 
Swedish  headform  is  shown  in  the  next  Fig.  22.  The  comparison  is 
very  good  except  at  a  velocity  of  30  feet  per  second  and  at  this  velocity 
the  difference  in  the  ct.  is  about  .  04.  It  is  also  seen  that  o.  increases 
with  increase  in  velocity;  this  trend  was  also  observed  for  the  case  of 
the  hemispherical  nose. 

Comparison  of  a.  with  -c0  :  Also  shown  on  Fig.  22  is  the  negative  value 
1  ^s 

of  the  pressure  coefficient  at  separation  which  is  obtained  with  the  help 

of  pressure  distribution  measurements  made  by  Rouse'  '  and  the 

observed  position  of  the  boundary  layer  separation  from  Fig.  15.  It 

should  be  noted  that  a.  is  closely  related  to  -c_  rather  than  -c  . 

i  Fg  pmin 

which  is  .  6  for  the  Swedish  headform.  The  difference  between  Oj  and 

-c  starts  with  a  value  of  about  0,075  at  a  velocity  of  26  feet  per 
ps 

second  and  gradually  decreases  to  zero  at  a  velocity  of  51  feet  per 


second. 
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6.  IMPORTANT  PHYSICAL  PARAMETERS  FOUND  IN  INCEPTION 

We  have  seen  on  both  test  bodies  inception  of  cavitation  to  be 
influenced  by  the  presence  of  the  laminar  separation.  However, 
inception  characteristics  were  different  in  the  sense  that  >  -  Cp^  for 
the  hemispherical  nose  and  <  -  Cp^  for  the  Swedish  headform  in  the 
same  Reynolds  number  range.  An  attempt  will  be  made  to  explain  the 
above  mentioned  difference  from  schlieren  observations  shown  in  Fig.  17. 

As  mentioned  earlier  in  section  4  one  of  the  principle  differences 
in  the  nature  of  real  fluid  flow  past  the  two  bodies  was  in  the  observa¬ 
tion  that  height  of  the  sep?~s>ted  layer  for  the  case  of  the  hemispherical 
nose  was  2  to  3  times  larger  than  for  the  case  of  the  Swedish  headform. 

This  difference  in  height  may  influence  the  magnitude  and 
frequency  of  the  pressure  fluctuations  in  the  transition  region  of  the 
separated  free  shear  layer  and  at  least  for  the  hemispherical  nose 
inception  was  seen  to  occur  in  this  region.  Gaster^0^  has  conducted 
some  qualitative  investigations  of  the  above  mentioned  point,  for  example 
he  observed  that  dominant  unstable  frequency  in  the  attached  boundary 
layer  upstream  of  separation  was  about  200  cycles  per  second,  but  on 
the  free  shear  layer  it  was  observed  to  be  about  1000  cycles  per  second. 
Detailed  measurements  of  pressure  fluctuation  needed  for  present 
analysis  are  not  available,  thus  the  reasoning  has  to  be  speculative 
in  nature. 

Present  observations  indicate  that  vertical  extent  of  the  laminar 
separated  region  significantly  controls  the  maximum  size  of  the  bubble 
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which  may  reside  in  the  region.  From  motion  pictures  of  the  electroly¬ 
sis  bubbles,  the  following  information  was  obtained  which  may  help 
understand  the  above  point: 


Model 

U 

m 

Tr 

Hm 

ins 

sec. 

in. 

Hem.  Nose 

30 

2.5X10-3 

6.5X10"3 

.012 

ft 

40 

2X10"3 

4.4X10"3 

.0084 

Swed.  Head. 

30 

1.1  xio“3 

6.6X10"3 

.005 

where  dm  is  the  maximum  diameter  of  the  bubble  stabilized  in  separated 
zone  in  inches,  Tr  is  the  residence  time  of  the  above  bubble  in  seconds, 
and  Hm  is  the  maximum  separated  height  from  Fig.  17. 

It  may  be  seen  that  dm  is  proportional  to  Hm.  It  should  be  seen 
that  for  vaporous  cavitation  to  occur  at  P  =  Pv  the  diameter  of  the  bubble 
should  be  greater  than  10"3  in.  Still  the  conditions  which  deter¬ 

mine  the  origination  of  the  microscopic  bubbles  without  electrolysis 
remains  to  be  found.  It  is  these  conditions  which  may  determine  the 
incipient  conditions  and  not  the  actual  magnitude  of  even  though 

latter  may  be  important  for  visual  method  of  determining  inception. 

From  Fig.  17,  for  both  models,  dm  decreases  with  increase  in  velocity 
and  thus  one  may  expect  to  decrease  with  increase  in  velocity,  but 
this  is  contrary  to  the  observations  of  Figs.  20  and  22.  From  this  one 
may  conclude  for  the  present  Reynolds  number  range,  that  the  height 
of  the  separated  layer  probably  plays  an  important  role  in  determining 
the  magnitude  of  pressure  fluctuation  in  the  transition  region  which  in 
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turn  might  determine  the  value  of  a..  But  this  may  not  hold  true  in 
generality. 

Under  any  circumstance,  separation  does  influence  the  process 
of  inception  significantly.  The  importance  of  c_  has  been  brought 
about  earlier  and  with  this  in  mind  a  new  cavitation  inception  coefficient 
is  hereby  proposed 

<cp,+°i> 

0*  =  --TS - r- 

1  (i-Cp  ) 
ys 

with  expected  additional  important  parameters  of  the  Reynolds  number 
based  on  momentum  thickness  at  separation  -Re?  ,  and  possibly  an 
additional  length  parameter,  say  H./8  where  H  is  a  relevant  length. 

It  is  tempting  to  think  that  this  length  might  be  proportional  to  maxi¬ 
mum  height  of  the  separated  layer  although  there  is  now  no  such  basis 
to  determine  such  a  length. 

Average  inception  data  from  present  studies  on  the  hemispherical 

/Q\ 

nose  (incipient  band  cavitation),  Swedish  headform,  Acosta's  data'  '  and 

/A)  $  ± 

Kermeen1 8  data'  1  are  plotted  in  terms  of  versus  Reg  in  Fig.  30. 

Kermeen' s  data  is  taken  from  Fig.  22  corresponds  to  the  values  when 

clear  steady  attached  cavity  was  first  formed  or  according  to  present 

terminology  when  band  type  inception  occurred.  Therefore,  it  may  be 

seen  now  that  when  known  state  of  cavitation  is  the  same  in  Kermeen' s 

and  the  present  case  the  corresponding  a.  values  agree  quite  well.  The 

method  of  obtaining  -cD  for  the  hemispherical  nose  in  Kermeen' s  case 

*s 

was  discussed  in  Section  5.1.  If  it  is  band  type  inception  and  the 
corresponding  value  of  a*  is  negative,  then  it  may  be  concluded  that  the 
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type  of  cavitation  is  gaseous  in  nature.  This  being  the  case  for  the 
hemispherical  nose  cavitation  observed  in  the  present  and  Kermeen's 
case. 

We  have  seen  that  cavitation  inception  process  is  influenced  by 

the  presence  of  a  laminar  separation  on  a  certain  class  of  bodies  and 

under  these  conditions  o.  may  be  considerably  different  from  -c  .  . 

1  p  min 

Furthermore,  it  is  also  observed  that  for  such  bodies  depending  on 

the  shape,  a.  may  be  closely  correlated  with  -c  .  But,  from  the 
1  ps 

present  experimental  program  no  explanation  is  found  for  the  observed 
size  effect  on  o.  for  geometrically  similar  bodies  of  certain  shape  for 
a  fixed  Reynolds  number.  Examples  of  such  dependence  are  well  known 
and,  in  particular.  Parkin  and  Holl^  have  observed  this  effect  on  a 
hemispherical  nose  and  a  1.5  calibre  ogive.  However,  it  is  believed 
that  proper  foundations  have  now  been  laid  by  the  present  experimental 
program  to  address  the  above  mentioned  size  effect  in  future  work. 
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7.  SUGGESTED  CORRELATION  TO  PREDICT  Cj  FOR 
BODIES  WITH  ATTACHED  BOUNDARY  LAYERS 


It  is  well  known  that  attached  laminar  boundary  layer  and  also 

laminar  separated  free  shear  layer  possess  a  point  of  instability  after 

which  the  outside  disturbances  grow  exponentially  in  amplitude.  The 

rate  of  growth  in  general  is  determined  by  the  local  Reynolds  number 

end  local  pressure  gradient  forces.  Such  spatial  growth  of  disturbances 

can  be  computed  for  an  attached  laminar  boundary  layer  with  the  help 

of  recent  computations  made  by  A.  M.  O.  Smith^\  It  has  been  pro- 
(41) 

posed  empirically  '  that  if  the  amplitude  of  the  initial  disturbances 

7  9 

grows  by  a  factor  of  exp  to  exp  ,  transition  takes  place  and  the  boundary 
layer  becomes  turbulent.  Detailed  calculations  with  the  help  of  results 
from  Ref.  (46)  in  determining  the  position  of  boundary  layer  transition 
using  above  criterion  is  explained  in  Appendix  B.  Using  such  calcula¬ 
tions  the  position  of  transition  has  been  estimated  as  a  function  o 1  Re^ 
for  a  1.5  calibre  ogive  with  cylindrical  afterbody.  Re^  crjt  for  1.5 
calibre  ogive  is  estimated  to  be  6.3X10^;  comments  regarding  the 
importance  of  Rej-j  crjt  have  been  made  earlier  in  section  4.  Thus, 
after  knowing  the  position  of  transition  and  pressure  distribution  on 
the  headform'  ,  one  can  determine  the  value  of  the  pressure  coeffi¬ 
cient  at  the  position  of  transition  (designated  by  c  ).  Variation  of 

”tr 

Cp^  with  Re^  for  15  calibre  ogive  is  shown  in  Fig.  B.5  of  Appendix  B. 

A  comparison  of  -  Cp^  and  desinent  cavitation  number  for  a 
two  inch  1.5  calibre  ogive  taken  from  Ref.  (5)  is  shown  in  Fig.  31.  These 
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results  show  the  close  correlation  of  o  ,  with  -  c-,  . 

a  Ptr 

The  above  finding  is  yet  to  be  verified  by  direct  experiments. 
However,  in  the  case  of  separated  flows,  inception  of  cavitation  was 
related  Lo  the  transition  region  of  the  separated  free  shear  layer,  a 
region  of  strong  pressure  fluctuations.  Thus,  the  comparison  afforded 
in  Fig.  31  suggests  that  the  somewhat  similar  fluctuations  in  the  transi¬ 
tion  region  of  the  boundary  layer  might  also  be  responsible  for  cavitation 
inception  on  bodies  not  having  a  laminar  boundary  layer  separation. 
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8.  DISCUSSION  OF  ELECTROLYSIS  STUDIES 


Swedish  Headform:  We  have  seen  in  Fig.  22  that  for  the  Swedish 

headform  was  considerably  lower  than  -  cpmj[n  and  in  some  cases  even 

lower  than  -  c  .  The  above  observation  indicates  that  just  prior  to 
Ps 

incipient  conditions  liquid  in  the  neighborhood  of  the  minimum  pressure 
point  is  under  an  actual  tension.  It  seemed  possible  when  this  work  was 
started  that  this  tension  could  not  be  supported  if  an  abundant  supply  of 
surface  nuclei  formed  by  electrolysis  were  to  be  made  available.  To 
investigate  this  aspect  of  cavitation  inception  such  nuclei  were  provided 
at  three  different  levels  of  tunnel  pressure  at  a  constant  velocity  of 
30  ft/ sec.  At  this  velocity  Oj  under  natural  conditions  was  found  to  be 
about  0.427,  and  photograph  (d)  of  Fig.  24  illustrates  this  condition. 

The  first  experiment  was  carried  out  at  a  cavitation  index  of 
a  =  0.852  which  is  greater  than  -  c-min>  “  cp  and  (Jj.  It  may  be  seen 
from  Fig.  24(a)  that  the  electrolysis  bubbles  do  not  show  any  growth, 
but  they  reside  in  the  separated  zone  much  longer  than  anywhere  else 

on  the  model.  From  motion  picture  studies  it  could  be  deduced  that 

-3  -4 

bubbles  larger  than  about  10  in.  in  diameter  stayed  less  than  10 

seconds  in  the  separated  zone  except  one  bubble  which  was  noticed  to 

be  entrained  in  the  reverse  flow  and  its  residence  time  was  found  to 

be  about  5  milliseconds. 

The  second  test  was  carried  out  for  a  value  of  a  -  0.544  which 
is  less  than  -  Cpmjn  but  greater  than  -  Cp  and  normal  from  (b) 
of  Fig.  24,  it  may  be  seen  that  electrolysis  bubbles  cavitated  readily 
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starting  from  the  minimum  pressure  point. 

The  third  and  last  value  a  =  .437  was  less  than  -  c  _ .  and  -  c 

pmin  ps 

but  is  still  greater  than  c^;  from  (c)  of  Fig.  24,  profuse  c  a  vita  don 
from  electrolysis  bubbles  starting  from  the  minimum  pressure  point  is 
apparent.  However,  the  presence  of  electrolysis  bubbles  did  not  trig¬ 
ger  a  steady  cavity  and  the  cavitation  observed  completely  disappeared 
as  soon  as  the  production  of  electrolysis  bubbles  was  stopped. 

Similar  results  were  observed  at  other  velocities  except  the 
physical  appearance  of  cavitation  from  electrolysis  bubbles,  as 
illustrated  in  the  photographs  of  Fig.  25,  changed  with  increase  in  velo¬ 
city.  In  Fig.  25(a),  at  30  feet  per  second  cavitation  occurred  as  isolated, 
large,  macroscopic  bubbles.  At  a  velocity  of  55  feet  per  second,  as 
shown  in  Fig.  25(c),  streaks  were  exhibited  and  at  an  intermediate 
velocity  of  40  feet  per  second  both  types  are  observed  simultaneously 
as  may  be  seen  in  Fig.  25(b). 

Hemispherical  Nose;  From  Fig.  23(a),  it  may  be  seen  that  for  CTXjj 
but  less  than  -  electrolysis  bubbles  do  cavitate  only  in  the 

transition  region  of  the  separated  free  shear  layer  as  was  the  case 
under  natural  incipient  conditions.  For  a  slightly  lower  o,  cavitation 
is  also  observed  in  the  minimum  pressure  area  as  shown  in  Fig.  23(b). 

In  Fig.  23(c)  a  photograph  of  cavitation  for  o<o^  under  natural  condi¬ 
tions  is  shown  and  it  may  be  observed  that  cavitation  is  in  the  separated 
layer  and  thus  begins  downstream  of  the  minimum  pressure  point  as 
previously  discussed.  The  next  photograph  Fig.  23(d)  shows  that 
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electrolysis  bubbles  readily  cavitate  starting  from  approximately  the 
minimum  pressure  point  at  the  same  value  of  o. 

Remarks:  It  may  be  said  from  the  above  studies  that  with  a  steady 

supply  of  sufficient  nuclei  one  can,  in  fact,  observe  steady  cavitation 

when  a  =  -  c  .  ,  and  this  relationship  may  be  used  to  determine 
pmin 

-  Cpm-n  for  complex  flows.  On  the  Swedish  headform,  just  prior  to 
inception,  a  considerable  portion  of  the  liquid  is  at  pressures  below 
vapor  pressure,  but  even  under  these  conditions  a  momentary  supply 
of  artificial  nuclei  did  not  trigger  a  steady  attached  cavity.  Only 
reduction  of  c  to  o.  produced  a  steady  attached  cavity. 
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9.  CONCLUSIONS  AND  SUMMARY 


A  schlieren  method  of  flow  visualization  widely  used  in  the  field 

of  compressible  aerodynamics  has  been  successfully  used  for  the  same 

purpose  in  water  tunnel  experiments.  This  method  has  shown  the 

hemispherical  nose  to  have  a  viscous  laminar  boundary  layer  separation 
5 

at  a  Rep  s  9.06  x  10  and  the  Swedish  headform  to  have  set  aration  at  a  Re^  = 
c 

8.78  X  10  .  Approximate  boundary  layer  stability  calculations  indicate 
ReDcrit"  5X106  for  the  hemispherical  nose  as  opposed  to  previous 
ReDcrit^  ^XIO^  indicated  by  Rouse^^  and  used  heretofor  by  all  re¬ 
searchers  in  the  area  of  cavitation  inception.  Similar  calculations 
indicate  Rej)crj^  1.5X10^  for  the  Swedish  headform.  The  length  of  the 
laminar  separated  zone  on  both  the  models  decreased  monotonically 
with  increase  in  Reynolds  number. 

On  the  hemispherical  nose,  cavitation  inception  in  the  form  of 
cavitation  bubbles  was  seen  to  first  occur  in  the  transition  region  of  the 
separated  free  shear  layer.  With  reduction  in  a  the  viscous  separated 
region  was  observed  to  fill  with  the  cavitation  bubbles  until  a  clear 
attached  cavity  was  formed.  Present  observations  regarding  cavitation 
inception  on  a  hemispherical  nose  strongly  indicate  that  similar  obser¬ 
vations  by  Parkin  were  within  the  viscous  separated  region  and  not 
submerged  in  the  attached  boundary  layer  as  was  formerly  thought. 

On  the  Swedish  headform  cavitation  inception  in  the  form  of  a 
clear  attached  cavity  was  seen  to  take  place  within  the  separated  region. 
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For  bodies  with  viscous  laminar  separation  the  important 

pressure  coefficient  to  correlate  c.  with  is  found  to  be  the  negative 

value  of  the  c  at  the  point  of  laminar  separation,  -cD  .  Other 
P  ps 

important  parameters  found  were  Reynolds  number  and  a  length  para¬ 
meter  which  is  a  measure  of  the  height  of  the  viscous  separated  layer. 

For  the  Swedish  headform  a  compared  favorably  with  -c_  and 

1  Fs 

inception  always  occurred  for  <7.£-cD  ;  the  equality  sign  to  hold  true 

at  high  Reynolds  numbers.  For  the  hemispherical  nose,  inception 

always  occurred  for  a.^-c  ;  it  is  believed  the  turbulent  fluctuations 

1  ps 

in  the  transition  region  of  the  separated  layer  and  gaseous  type  of  cavi¬ 
tation  may  be  responsible  for  the  above  observation. 


From  present  experiments  no  explanation  could  be  offered  for 
the  observed  size  effect  on  a.  for  geometrically  similar  bodies  at  fixed 
Reynolds  numbers.  The  schlieren  method  may  be  useful  in  shedding 
light  on  this  question  in  future  work. 


Computation  of  the  pressure  coefficient  at  transition,  c  ,  on 

ptr 

1.5  calibre  ogive  for  Rep>Re^  cr-t  when  viscous  laminar  separation 


is  absent,  has  shown  that  a,  i  -c  . 

d  ptr 

expe  rimental  ly . 


This  finding  has  yet  to  be  verified 


LIST  OF  FIGURE  CAPTIONS  -  PART  I 


Fig.  1.  Cavitation  inception  number  versus  velocity  at  maximum 
acoustic  pressure  for  a  series  of  6  hemispheres  with 
cylindrical  afterbody  of  diameter-D  (Ref.  4). 

Fig.  2.  Results  of  a  comparative  cavitation  inception  tests  on  a 

modified  ellipsoidal  headform  sponsored  by  the  International 
Towing  Tank  Conference.  Taken  from  Fig.  8  of  Ref.  (7). 

Fig.  3.  Schematic  drawing  showing  the  schlieren  set-up  in  the  Free 
Surface  Water  Tunnel. 

Fig.  4.  Schematic  diagram  showing  the  mounting  of  a  two  inch  hemi¬ 
spherical  nose  in  the  Free  Surface  Water  Tunnel.  Additionally 
the  different  places  of  temperature  measurements  to  obtain 
heat  transfer  data  are  shown. 

Fig.  5.  Schematic  diagram  of  the  schlieren  set-up  in  the  High  S’-e-jd 
Water  Tunnel. 

Fig.  6.  Illustration  of  light  source  side  of  the  schlieren  system.  On 
the  far  right  is  the  spark  gap  source,  followed  by  deflecting 
mirror  and  columnating  lens.  The  white  box  on  the  far  left 
is  the  steady  light  source.  A  piano  convex  correction  lens  is 
shown  taped  on  the  tunnel  window. 

Fig.  7.  Illustration  of  camera  side  of  the  schlieren  system.  In  this 
figure  the  first  lens  is  the  focusing  lens,  followed  by  the  cut¬ 
off  knife  edge  placed  just  before  the  camera  box  itself.  The 
camera  box  is  80  in.  long. 

Fig.  8.  A  drawing  showing  the  details  of  placing  the  electric  heater 
inside  the  Swedish  headform.  The  dark  shaded  area  is  the 
heater. 

Fig.  9.  Diagram  showing  the  geometrical  details  of  the  Swedish 
headform. 


Fig.  10. 
Fig. 11. 

Fig.  12. 

Fig.  13. 


Fig. 14. 

Fig.  15. 
Fig.  16. 
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Diagram  showing  the  geometrical  details  of  the  hemispherical 
nose. 

A  schlieren  photograph  showing  the  laminar  separation  bubble 
on  a  two  inch  hemispherical  nose  at  a  Reynolds  number  of 
105.  The  flow  is  from  right  to  left.  Also  noticed  is  the 
gradual  break-up  of  the  separated  free  shear  layer  and  its 
subsequent  reattachment  to  the  model  surface. 

Photographs  illustrating,  that  the  position  or  extent  of  laminar 
separation  is  not  influenced  by  the  change  of  model  surface 
temperature.  Top  photograph  is  with  surface  temperature 
about  10  °F  below  the  free  stream  water  temperature  and  the 
bottom  is  with  surface  temperature  about  10°F  above.  Both 
photographs  are  of  a  two  inch  hemispherical  nose  at  Re  =  2.5X10 

Schlieren  photographs  illustrating  the  change  in  extent  of 
laminar  separation  bubble  with  increase  in  tunnel  velocity  or 
Reynolds  number  for  both  the  test  models.  The  top  two 
photographs  (a)  and  (g)  are  taken  with  no  model  heating  at  a 
tunnel  velocity  of  40  ft/sec.  (a)  Re^  =  5.85  X  10^, 

(b)  ReD=  2.93X  105,  (c)  ReD  =  4.  39X  105,  (d)  ReD=  5.85 X  105, 

(e)  ReD=7.3lX  105,  (f )  ReD=  8.78X  105,  (g)  ReD=  6.04X  10S, 

(h)  ReD=3.02xi05,  (i)  ReD=4.53X105,  (j)  ReD  =  6.04xl05 
(k)  ReD=7.55  X  105,  (f )  ReD  =  9. 06  X  105. 

Variation  of  laminar  boundary  layer  separation  position  on 
the  hemispherical  nose  with  tunnel  velocity.  Also  shown  is 
the  effect  of  temperature  change  on  the  position  of  separation. 

Variation  of  the  position  of  boundary  layer  separation  on  the 
Swedish  headform  with  change  in  tunnel  velocity. 

Variation  of  the  length  of  the  laminar  separated  bubble  estima¬ 
ted  from  schlieren  photographs  with  tunnel  velocity  for  both 
the  test  models. 


I 


I 


Fig.  17. 


Fig.  18. 


Fig.  19. 


Fig.  20. 


Fig.  21. 


Fig.  22. 


Fig.  23. 


Fig.  24. 
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Variation  of  the  maximum  height  of  the  laminar  separated 
bubble  indicated  by  the  schlieren  photographs  with  velocity 
for  both  the  test  models. 


Sequence  showing  the  progressive  development  of  cavitation 
on  the  hemispherical  nose  in  the  laminar  boundary  layer 
separated  region  at  a  tunnel  velocity  of  40  ft/ sec.  corres¬ 
ponding  Re^  =  6.04X  10^.  The  dark  patches  above  the  model 
surface  are  cavitating  areas. 

Partial  sequence  showing  the  progressive  development  of 
cavitation  on  the  hemispherical  nose  in  the  separated  region 

5 

at  a  tunnel  velocity  of  60  ft/sec.  corresponding  Re^  =  9.06  x  10 

Variation  of  cavitation  inception  number  wtih  tunnel  velocity 
on  the  hemispherical  nose  for  two  types  of  incipient  conditions. 
Also  included  are  Kermeen's  data. 


Variation  of  cavitation  inception  number  for  the  Swedish 
headform  with  tunnel  velocity.  Inception  was  always  of  the 
"band"  type. 

Comparison  of  present  average  cavitation  inception  data  with 

(9) 

that  obtained  by  Acosta'  '  on  the  same  Swedish  headform  and 
also  with  the  negative  value  of  the  pressure  coefficient  at 
laminar  separation. 


Sequential  photographs  illustrating  the  effect  of  introducing 
artificial  surface  nuclei  by  electrolysis  on  the  hemispherical 
nose  at  different  tunnel  pressures.  The  photograph  (c)  is 
under  natural  cavitating  conditions.  U  =  28  ft/sec. ,  a  =  .72, 

"cpmin=*74,  "CPS~*63‘ 

Photographs  illustrating  the  effect  of  introducing  artificial 
surface  nuclei  by  electrolysis  on  the  Swedish  headform  at 
different  tunnel  pressures.  Cavitation  observed  in  photo¬ 
graphs  (b)  and  (c)  would  not  be  present  at  the  corresponding 
cavitation  numbers  without  electrolysis.  Photograph  (d) 


r 


Fig.  25. 

Fig.  26. 

Fig.  27. 

Fig.  28. 

Fig.  29. 

Fig.  30. 
Fig.  31. 
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shows  the  natural  incipient  cavitation.  U=  31  ft/sec. ,  Cj=  .4, 
-cpmin=-6'  -<V'46- 

Photographs  illustrating  the  change  in  the  physical  appearance 
of  cavitation  produced  due  to  electrolysis  bubbles  with  increase 
in  tunnel  velocity  on  the  Swedish  headform.  For  comparison 
physical  appearance  of  cavitation  under  natural  conditions  is 
shown  in  the  photograph  (d). 

Sequential  photographs  showing  two  different  types  of  up¬ 
stream  motion  of  the  microscopic  bubbles  within  the  separated 
layer  on  the  hemispherical  nose.  The  free  stream  flow  is  from 
left  to  right.  The  dark  patches  at  the  downstream  end  of  the 
photographs  is  the  macroscopic  cavitation  readily  observable 
to  the  naked  eye.  U  =  40  ft/sec. ,  o  =  .74. 

Sequential  photographs  illustrating  the  nature  of  desinent 
cavitation  on  the  Swedish  headform.  Last  traces  of  cavitation 
are  seen  to  exist  at  the  end  of  the  separated  layer.  U  =  40 
ft/ sec. ,  a=  .46.  The  top  schlieren  photograph  is  at  U=  30 
ft/  sec.  The  free  stream  flow  is  from  left  to  right. 

Variation  of  the  length  of  the  laminar  separated  bubble  for 
the  two  test  models  in  non-dimensional  form  with  Reynolds 
number.  Also  shown  is  similar  data  obtained  by  Gaster^^^ 
in  wind  tunnel  studies. 

Variation  of  non-dimensional  maximum  height  of  the  laminar 
separated  bubble  indicated  by  the  schlieren  photographs  for 
both  the  test  models  with  Reynolds  number. 

Variation  of  cavitation  inception  coefficient  introduced  in 
present  work  with  Reynolds  number  for  different  models. 

Comparison  of  incipient  cavitation  index  for  1.5  calibre 
ogive  with  the  negative  value  of  the  pressure  coefficeint  at 
the  predicted  point  of  transition  for  Rep>Rep  cr^» 


TUNNEL  VELOCITY -U  IN  FEET  PER  SECOND 

Fig.l.  Cavitation  inception  number  versus  velocity  at  maximum 
acoustic  pressure  for  a  series  of  6  hemispheres  with  cylindrical 
afterbody  of  diamcter-D  (Ref. 4.  ). 


COMPARATIVE  INCEPTION  STUDIES  ON 
A  MODIFOID  ELLIPSOIDAL  HEADFORM 
TAKEN  FROM  FIG.  8  OF  REF.  (7) 
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,  5.  Schematic  diagram  of  theschlieren  set-up  in  the  High  Speed 
er  Tunnel. 


schlieren  system.  In  this  figure,  the  first 
lens  is  the  focusing  lens,  followed  by  the 
cut  off  knife  edge  placed  just  before  the 
camera  box  itself.  The  camera  box  is  80  in. 
long. 
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MATERIAL -STAINLESS  STEEL 
SURFACE  FINISH -SMOOTH 
AXISYMMETRIC 


Fig. 9.  Diagram  showing  the  geometrical  details  of  the  Swedish 
headform. 


MATERIAL-  BRASS 
SURFACE  FINISH- SMOOTH 
AXISYMMETRIC 

Fig.  10.  Diagram  showing  the  geometrical  details  of  the  hemispheri¬ 
cal  nose. 


SWEDISH  HEADFORM  -  D  =  1,755  in.  HEMISPHERICAL  NOSE  -  D=  1 . 813  in 
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Fig.  17.  Variation  of  the  maximum  height  of  the  laminar  separated 
bubble  indicated  by  the  schlieren  photographs  with  velocity  for  both 
the  test  models. 


(a)  a  =.772 


(b)  a  =.70 


— 


MODEL  IN  SILHO 


gEHEEl 


Fig.  18.  Sequence  showing  the  progressive  development  of  cavitation 
on  the  hemispherical  nose  in  the  laminar  boundary  layer  separated 
region  at  a  tunnel  velocity  of  40  ft/  sec.  corresponding  Rep  =  6.04  x  105 
The  dark  patches  above  the  model  surface  are  cavitating  areas. 


\ 


(d)  o=  .85 


(e)  CT=  .732 

Fig.  19.  Partial  sequence  showing  the  progressive  development  of 
cavitation  on  the  hemispherical  nose  in  the  separated  region  at  a 
tunnel  velocity  of  60  ft/sec.  corresponding  Re^  9.06  x  10*. 
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TUNNEL  VELOCITY- U  IN  FEET  PER  SECOND 

Fig.  20.  Variation  of  cavitation  inception  number  with  tunnel  velocity- 
on  the  hemispherical  nose  for  two  types  of  incipient  conditions.  Also 
included  are  Kermeen's  data. 
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(a)  a =.731 


(b)  a =.722 


k 


5  ?VJCAVITy  :_  '  . 

NATURALCAVITATI01^^^^^^^^^^^H(c)ct=^597 


(d)  o  =  .597 

Fig,  23.  Sequential  photographs  illustrating  the  effect  of  introducing 
artificial  surface  nuclei  by  electrolysis  on  the  hemispherical  nose  at 
different  tunnel  pressures.  The  photograph  (c)  is  under  natural  cavi- 
tating  conditions.  U=28ft/sec. ,  q=,72,  -c  ^^=.74,  -c  i  .63 


Fig.  24.  Photographs  illustrating  the  effect  of  introducing  artificial 

surface  nuclei  by  electrolysis  on  the  Swedish  headform  at  different 

tunnel  pressures.  Cavitation  observed  in  photographs  (b)  and  (c)  would 

not  be  present  at  the  corresponding  cavitation  numbers  without  electroly 

sj.s.  Photograph  (d)  shows  the  natural  incipient  cavitation.  U=  31  ft/ sec. 

■Cfe  =  .4.  -  c  .  =  .6,  -  c  =,46. 

*  prrnn  '  ps 
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Fig.  25.  Photographs  illustrating  the  change  in  the  physical  appearance 
of  cavitation  produced  due  to  electrolysis  bubbles  with  increase  in 
tunnel  velocity  on  the  Swedish  headform.  For  comparison  physical 
appearance  of  cavitation  under  natural  conditions  is  shown  in  the 
photograph  (d). 
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'Av\y, 


0.00225 


g.  26a 


Fig.  26b 


Fig.  26.  Sequential  photographs  showing  two  different 
types  of  up  stream  motion  of  the  microscopic  bubbles 
within  the  separated  layer  on  the  hemispherical  nose. 
The  free  stream  flow  is  from  left  to  right.  The  dark 
patches  at  the  downstream  end  of  the  photographs  is  the 
macroscopic  cavitation  readily  observable  to  the  naked 
eve.  U=  40  ft/sec. .  a  =.74. 


0.00225 


Fig.  27.  Sequential  photographs  illustrating  the  nature  of  desinent 
cavitation  on  the  Swedish  headform.  Last  traces  of  cavitation  are  seen 
to  exist  at  the  end  of  the  separated  layer.  U=  40  ft/sec. ,  a=  .46.  The 
top  schlieren  photograph  is  at  U=  30  ft/sec.  The  free  stream  flow  is 
from  left  to  right. 
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HEMISPHERICAL  NOSE-  D  =  1.813  INCHES 


oc  of 

u  •> 
li->-  • 

o<=> 

-J  I 
Z  2 

Q£o 

Wq< 

OzK 

0-D< 

ffiUJ 
“  </) 


£.*<6 
tz  a> 

U  QC 

3  - 

fO  UJ  Z 

>  2 

z  5 
o  o: 
< 
O  D- 

8  y 

0(/5  ^ 

o  <r 

ft  to 

g9  LU 

1  z 
3  * 

oZ2 
o  X 


v  c 

g*z 

u  c  « 


cd  «  p 

ES* 

g_.  o 

OT3  d 
•h  w  r 

CD  4-» 

d  «  <D 
8 

111 

•  «  ? 
RH  . 
O  ^  ” 

fl*2  « 

<m 

o  o 
fill  fi 

o  «  M 
H  4J  4J 

d  « 
flj  M  « 

TJ  4  ** 

S  8* « 

.  h  A 

Sg* 

•  •S'0 

60  g  u 
•H  3  JJ 


8^/X0U,H 

•N0llVdVd3S  IV  SS3N»0IHi  WfUN3W0W  01 

3iaana  asivavdss  do  ihogh  wnwixvw 


91- 


00 

CD 

*■ 

CM 

o 

CM 

sr 

CD 

o 

o. 

q 

q 

q 

q 

q 

• 

o 

o 

d 

o 

o 

o 

i 

o 

i 

-lN3IDIdJ300  NOIld30NI  N0I1V1IAV0 


PART  II 


A  SEMI-  EMPIRICAL  METHOD  TO  PREDICT 
CAVITATION  SEPARATION  ON  SMOOTH  BODIES 
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1.  INTRODUCTION 


In  Part  I,  we  have  considered  the  problem  of  cavitation  inception 
on  bodies  which  possess  laminar  boundary  layer  separation.  In  the 
present  section  we  are  concerned  with  the  problem  of  subsequent  cavi¬ 
tation  development  on  such  bodies  with  the  reduction  of  a  sufficiently 
below  that  of  a..  In  particular,  we  are  interested  in  viscous  effects  on 
the  beginning  position  of  the  developed  cavity,  which  we  may  call  the 
position  of  cavitation  separation. 

In  order  to  predict  analytically  the  force  coefficients  on  bodies 
with  fully  developed  cavities,  a  knowledge  of  the  position  of  cavitation 
separation  is  essential.  Bodies  with  sharp  corners  such  as  discs, 
wedges  and  sharp-edged  hydrofoils  at  large  angles  of  attack  possess 
cavitation  separation  whose  position  is  known  apriori  and  force  coeffi¬ 
cients  for  such  bodies  can  be  predicted  quite  well  within  the  potential 

(47) 

flow  approximations  But,  for  the  case  of  smooth  bodies  the  condi¬ 

tion  of  "smooth  separation"  is  used  to  predict  the  position  of  cavitation 
detachment'  .  The  above  mentioned  condition  does  not  take  into 
account  the  possibilities  of  real  fluid  effects  on  the  position  of  cavitation 
separation.  Recent  experimental  studies^®’  ^  have  shown  that  the 
position  of  cavitation  separation  on  smooth  bodies  does  depend  on 
Reynolds  number  or  is  influenced  by  viscous  effects.  Further,  it  was 
also  observed  that  the  experimentally  found  position  of  cavitation 
separation  lies  considerably  downstream  than  that  predicted  by  the 


condition  of  smooth  separation.  On  a  sphere  the  difference  was  of 

4 

the  order  of  40°  at  an  Re^  and,  the  difference  decreased  with 

increase  in  Re^. 

Similar  results  were  noted  for  the  case  of  a  bi-convex  hydro¬ 
foil^.  On  this  hydrofoil,  linearized  theory  with  cavitation  separation 
positions  determined  by  smooth  separation  condition  predicted  positive 
lift  coefficients  at  small  positive  angles  of  attack.  Experiments,  on 
the  contrary,  showed  negative  lift  coefficients  at  the  same  angles  of 
attack.  This  discrepancy  was  directly  related  to  the  erroneous  pre¬ 
diction  of  cavitation  separation  position  by  the  smooth  separation 
criterion.  This  example  illustrates  the  extremely  important  role  the 
position  of  the  point  of  cavitation  separation  has  in  determining  overall 
forces,  yet  very  little  attention  has  been  devoted  to  this  subject. 

In  the  latter  study^  it  is  mentioned  that  two  different  types  of 
cavitation  separation  are  observed  on  a  smooth  body,  namely  a)  vis¬ 
cous  laminar  cavitation  separation  and  b)  nucleate  cavitation  separation. 
Physical  appearance  of  these  two  types  of  reparation  is  very  different 
as  illustrated  by  the  two  photographs  of  Fig.  1.  Brennen^®^  has  observed 
the  type  (a)  separation  on  a  sphere  even  at  Reynolds  number  of  10^. 

From  present  studies  of  Part  1,  we  have  a  better  idea  of  what  is  physi¬ 
cally  involved  in  the  two  types  of  separation.  It  must  be  concluded 
that  viscous  laminar  cavitation  separation  is  observed  only  on  bodies 
which  possess  a  laminar  boundary  layer  separation  under  non-cavitating 
conditions.  For  such  bodies  it  was  observed  in  Part  I  that  inception  of 
cavitation  does  take  place  first  within  the  separated  layer  and  if  for 
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some  reason  it  did  not,  as  was  the  case  when  electrolysis  bubbles 
were  introduced,  nucleate  type  separation  was  observed  on  the  same 
body  and  at  the  same  Reynolds  number  as  shown  by  Fig.  24(c)  and 
Fig.  24(d)  of  Part  I,  Thus,  it  is  apparent  that  the  type  of  cavitation 
separation  to  be  expected  on  a  smooth  body  is  dependent  on  the  com¬ 
plete  history  of  cavitation  development. 

Part  II  is  iniended  to  take  a  deeper  look  at  the  nature  of  real 
fluid  flow  in  the  neighborhood  of  viscous  laminar  cavitation  separation 
with  the  help  of  the  schlieren  method  as  flow  visualization  technique. 
Furthermore,  an  attempt  will  be  made  to  develop  a  semi-empirical 
method  of  predicting  cavitation  separation  on  bodies  which  possess 
laminar  boundary  layer  separation  under  fully  wetted  conditions  with 
the  hope  that  the  present  method  will  be  superior  to  that  of  using  the 
smooth  separation  condition. 
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2.  EXPERIMENTAL  PROCEDURE 

The  same  experimental  methods  described  in  section  2  of  Part 
I  were  used  for  present  work.  Schlieren  silhouette  photographs  of  the 
fully  developed  cavity  were  taken  with  the  spark  gap  light  source  for 
both  the  models,  at  various  tunnel  velocities  and  cavitation  numbers. 
Such  a  photograph  is  illustrated  in  Fig.  2  for  the  Swedish  headform 
taken  at  a  tunnel  velocity  of  25.6  ft/sec.  and  c  =  0.39.  It  should  be 
noted  from  this  photograph  that  the  cavity  separation  is  preceded  by 
the  laminar  boundary  layer  separation.  The  streamwise  distance  be¬ 
tween  the  two  separations  is  denoted  by  "X".  With  the  help  of  such 
photographs  it  was  possible  to  determine  the  position  of  the  boundary 
layer  separation,  the  position  of  cavitation  separation  and  hence  the 
distance  between  the  two,  X. 
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3.  DISCUSSION  OF  EXPERIMENTAL  RESULTS 

Variation  of  the  angular  position  of  the  cavity  separation  with 
tunnel  velocity  and  cavitation  number  for  the  hemispherical  nose  is 
shown  in  Fig.  3.  From  this  Fig.  3  it  may  be  seen  that  the  position  of 
the  cavity  separation  shifts  upstream  by  about  12°  with  a  decrease  of 
o  from  .65  to  .25,  and  it  shifts  downstream  by  about  8°  for  a  decrease 
in  tunnel  velocity  from  55  ft/sec.  to  25  ft/ sec.  Similar  results  have 
been  noted  for  a  sphere^)  jn  the  same  Reynolds  number  range.  Fig¬ 
ure  4  shows  the  variation  of  the  angular  position  of  the  boundary  layer 
separation  which  precedes  the  cavity  separation  on  the  hemispherical 
nose  as  a  function  of  tunnel  velocity  and  a.  It  is  important  to  note  from 
this  Fig.  4,  th*t  unlike  the  position  of  the  cavity  separation,  the  position 
of  the  boundary  layer  separation  preceding  it  is  independent  of  the 
tunnel  velocity. 

Variation  of  the  streamwise  position  of  the  cavity  separation 
and  the  boundary  layer  separation  preceding  it  for  the  Swedish  headform 
as  function  of  the  tunnel  velocity  and  o  is  shown  in  Figs.  5  and  6.  It 
may  be  noticed  that  the  results  are  similar  to  that  observed  for  the 
hemispherical  nose  and  discussed  earlier.  From  the  way  results  are 
presented  in  Ref.  (9)  it  is  not  possible  to  make  a  direct  comparison  be¬ 
tween  the  presently  observed  positions  of  the  ca\ity  separation  and  that 
/o\ 

obtained  by  Acosta'  \  But,  from  Fig.  5  it  may  be  seen  that  the  general 
agreement  between  the  two  results  mentioned  above  is  fairly  good. 
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The  ratio  of  X  to  the  diameter  of  the  model,  D,  is  plotted  against 
a  for  various  tunnel  velocities  for  the  hemispherical  nose  in  Fig.  7  and 
for  the  Swedish  headform  in  Fig.  8.  From  these  figures  it  may  be  noted 
that  X/D  increases  with  decrease  in  tunnel  velocity,  but  is  almost  con¬ 
stant  with  change  in  a  or  at  least  did  not  show  any  consistent  depen¬ 
dence  with  a  for  a  fixed  velocity.  To  elucidate  the  latter  point  more 
clearly,  the  variation  of  X/D  is  plotted  against  o  in  Figs.  9  and  10  for 
both  the  models,  but  at  only  two  velocities  for  which  the  maximum 
amount  of  data  was  obtained.  These  two  figures  illustrate  the  point 
that  X/D  is  almost  independent  of  a,  even  though  for  the  Swedish  head- 
form  at  a  tunnel  velocity  of  40  ft/ sec.  the  scatter  is  considerable.  It 
should  be  noted  that  in  general  the  values  of  X/D  were  quite  small;  the 
maximum  being  0.07  and  the  minimum  of  0.01  thus,  some  scatter  in 
the  measurement  of  these  small  distances  should  be  expected. 

With  the  earlier  observation  of  the  independence  of  X/D  with  o, 

it  is  justifiable  to  take  an  average  value  of  X/D  for  all  the  meaningful 

$ 

cavitation  numbers  and  this  average  will  be  denoted  by  X  /D.  Varia- 
tion  of  X  /D  with  tunnel  velocity  for  both  the  models  is  shown  in  Fig.  1 1. 
Each  data  point  in  Fig.  11  represents  the  averaged  value  of  X/D  along 
the  horizontal  axis  of  the  Figs.  7  and  8  at  the  corresponding  velocity. 

It  is  evident  that  one  can  determine  the  position  of  cavitation 
separation  once  the  position  of  the  boundary  layer  separation  preceding 
it  and  the  quantity  X  is  known.  On  present  models,  observations 
strongly  suggest  that  the  position  of  the  boundary  layer  separation 
preceding  the  cavity  separation  is  dependent  only  on  a,  but  is 
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independent  of  the  tunnel  velocity  and  the  quantity  X  is  only  dependent 
on  the  tunnel  velocity,  but  is  independent  of  o.  Next  an  attempt  will  he 
made  to  generalize  the  above  mentioned  observations  based  on  physical 
arguments  to  obtain  a  correlation  using  the  present  results.  Then  this 
correlation  will  be  applied  to  predict  cavitation  separation  on  a  sphere 
and  a  cylinder. 


i 


4.  A  CORRELATION  FOR  BOUNDARY  LAYER 
SEPARATION  WITH  CAVITATION 


In  Part  I  we  have  seen  the  importance  of  the  pressure 
coefficient  at  boundary  layer  separation  in  cavitation  inception;  it  was 
also  observed  that  the  beginning  of  the  macroscopic  cavitation  com¬ 
menced  only  for  cavitation  numbers  lower  than  the  negative  value  of 
the  pressure  coefficient  mentioned  above.  Since  fully  developed 
cavitation  occurs  for  lower  cavitation  numbers  than  for  the  incipient 
conditions  it  was  found  appropriate  to  use  the  difference  of  the  values 
as  the  significant  parameter  which  may  be  expressed  in  a  normalized 
form  as  [l+c/(Cps)N  q  3.  This  term  will  be  denoted  by  a',  and  it  may 
be  thought  of  as  a  measure  of  the  extent  of  cavitation. 

As  may  be  seen  from  Fig.  18(c)  of  Part  I,  the  presence  of  the 
macroscopic  cavity  shifts  the  location  of  the  original  boundary  layer 
separation  upstream  to  a  new  position.  It  seems  appropriate  that  one 
should  consider  only  this  shift  as  a  physically  important  quantity.  Fur¬ 
ther  it  was  observed  that  for  the  same  change  in  the  parameter  o',  the 
above  mentioned  shift  was  two  to  three  times  larger  on  the  Swedish 
headform  than  the  hemispherical  nose.  It  is  believed  that  this  is  due  to 
the  fact  that  the  hemispherical  nose  possesses  a  steeper  pressure 
gradient  approaching  the  point  of  boundary  layer  separation  (seee.g. ,  Fig. 
A.  1,  Appendix  A).  To  take  into  account  the  above  observations  it 
seemed  appropriate  to  normalize  the  shift  in  the  boundary  layer 
separation  by  the  distance  between  the  minimum  pressure  point  and 
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the  position  of  the  boundary  layer  separation.  The  latter  distance  was 
taken  to  be  that  existing  under  non-cavitating  conditions.  Symbolically, 
the  non-dimensional  groups  found  were, 


cy/=  [1  +  CT/(CP8)N.  c.  ]  and  i=  JZZ 


B.L.S./D)N.C.  "  ^B.L.s/^g 
B.L.S./D)N.C.  (scpmin/D)N.C.. 


A  plot  of  o' versus  $  for  the  present  models  is  shown  in  Fig.  12  and  it 
may  be  said  that  the  data  collapses  fairly  well.  In  the  construction  of 
Fig.  12  the  following  information  was  utilized: 

a)  (8tj  r  e  /D)«t  ~  =0.83  for  the  Swedish  headform  and  is 

In,  L, 

obtained  from  Fig.  15  of  Part  I.  Corresponding 

(cp  )N  c  =  -.46  from  Fig.  A.  1  of  Appendix  A. 

8 

b)  (Sg  L  g  /D)^  ^  =0.758  for  the  hemispherical  nose  and  is 
obtained  from  Fig.  14  of  Part  I.  Corresponding 

(cp  )N  c  =  -.63  from  Fig.  A.l  of  Appendix  A. 

c)  <8B.L.S./D|N.C.'(Bcpmin/D,N.C.  *0.2  for  the  Swedish 
headform  and  is  obtained  from  (sc  T  e  /D).t  ~  ±  0.83 

and  Fig.  A.l  of  Appendix  A  where  *8cpmli/D,N.C.  ~0-83‘ 

d>  (8B.L.S.  /“'Re.  *  <8  cpmin/D)N.C.-0-078  SoT  ,he 

hemispherical  nose  and  is  obtained  from  (Sg  ^  g  /D)^  ^  ±  0.758 

and  Fig.  A.l  of  Appendix  A  where  (s  cpmta/D>N.C.i0-68- 
A  straight  line  approximation  for  the  data  shown  in  Fig.  12  gives 
♦  =  2.37[l  +  a/(cp  )K  C  ]  (1) 

It  is  clear  that  more  data  is  needed  to  extrapolate  the  present  data 
with  a  straight  line.  But,  later  we  will  see  that  when  this  straight 


line  approximation  is  used  to  compute  the  position  of  cavitation  separa¬ 
tion  on  a  sphere  at  o'  as  high  as  .9  good  agreement  is  obtained  with 
experimentally  observed  values. 

Therefore,  for  given  information  under  non-cavitation  conditions 
and  a,  it  is  possible  to  determine  the  position  of  the  boundary  layer 
separation  which  precedes  cavitation  separation  with  the  help  of 
Equation  (1). 

I 

I 

( 
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5.  A  CORRELATION  FOR  THE  DISTANCE 
BETWEEN  THE  TWO  SEPARATIONS 


The  physical  appearance  of  cavitation  occurring  on  present 

models  as  represented  by  Fig.  1(b)  has  striking  similarity  with  the 

physical  appearance  of  cavitation  (in  the  form  of  air  entrainment) 

occurring  in  the  diverging  sections  of  oil  film  lubrication  bearings  as 

(49) 

illustrated  by  Taylor  ,  Even  though  the  two  flow  configurations  are 
quite  different  especially  in  terms  of  inertial  effects,  the  similarity  is 
provided  from  the  fact  that  cavitation  on  the  present  models  occurs  in 
the  separated  zone  of  the  viscous  laminar  boundary  layer  (see  Fig.  18 
of  Part  I)  where  inertial  effects  may  be  quite  small  as  compared  to  the 
inertial  effects  in  the  outer  flow.  Pearson^^  has  conducted  some 

(49) 

experiments  which  are  similar  in  nature  to  that  conducted  by  Taylor'  , 
but  with  a  simpler  apparatus.  The  apparatus  consisted  of  a  symmetric 
wedge  with  divergence  angle  of  1/20  placed  over  a  flat  plate  with  a 
minimum  clearance  of  0.004  in.  and  the  region  between  the  two  was 
filled  with  different  types  of  liquids  possessing  high  viscosities  com¬ 
pared  to  water.  The  symmetric  wedge  was  terminated  at  x=  +lin.  and 
x  =  -Hn.  in  the  notation  of  Fig.  13.  It  was  drawn  at  measured  constant 
speed  perpendicular  to  planes  x=  const.  This  resulted  in  the  suction 
of  an  air  cavity  between  the  diverging  section  of  the  apparatus  and  the 
position  of  the  meniscus  was  determined  with  flash  photographs.  The 
air  cavity  possessed  a  characteristic  instability  in  the  z-direction  with 
alternate  presence  of  air  fingers  and  fully  wetted  regions.  This  is  very 
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similar  to  the  physical  appearance  of  cavitation  on  present  bodies  as 
shown  in  Fig.  1(b). 

For  the  problem  studied  by  Pearson,  the  important  dimension- 

|UU0 

less  group  was  found  to  be  -y  with  h^  as  the  pertinent  length  para- 
|iUn 

meter,  y  is  commonly  known  as  Taylor -Saffman  parameter  and  is 

the  ratio  of  Weber  number  to  Reynolds  number,  both  based  on  Uq.  The 

physical  conditions  under  which  cavitation  occurs  on  the  present  models 

is  shown  in  Fig.  14,  and  the  similarity  with  conditions  shown  in  Fig.  13 

is  very  suggestive.  With  this  in  mind  the  important  dimensionless 

group  was  taken  to  be  Taylor-Saffman  parameter  with  Weber  number 

based  on  local  velocity  at  separation  and  Reynolds  number  based  on 

free  stream  velocity.  In  the  present  symbolic  notation,  the  parameter 

would  be  [1  -  ( c  )N  £  ].  The  length  parameter  equivalent  to  h^ 

*  * 

was  chosen  to  be  the  momentum  thickness  at  separation. 

As  discussed  earlier  it  was  noted  that  X  was  independent  of  c 
and  this  allows  us  to  choose  the  pertinent  variables  in  determining  > 
at  any  cavitation  number  and  the  most  useful  condition  to  choose  would 


be  that  under  non- cavita ting  conditions  or  for  C>>1.  Thus,  ths  non- 

jgg  Li  U 

dimensional  groups  chosen  were  X  /(0fl)N  q  and-y  [l-(cp  )N  c  1  A 
plot  of  these  quantities  for  the  present  models  is  shown  in  Fig.  15.  The 


value  of  (0s)jj  q  was  estimated  by  Thwaite's  approximate  method  as 
outlined  in  Appendix  A.  Also  included  in  Fig.  14  are  results  obtained 


by  Pearson  (see  Fig,  11  of  Ref.  (50) );  b/hg  in  Pearson's  experiments 
is  equivalent  to  X* /(0g)N  q  for  the  present  experiments  and  iiHll  is 
equivalent  to  [l-  (cp  Jjq.c.  ^  **  maY  b®  note<*  from  Fig.  15  that  the 
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agreement  between  prrnrnf remits  Tnfl-ttnt-nf  Pearson’s  may  be  con- 
sidered  quite  good  especially  for  low  values  of  the  Taylor -Saffman 
parameter.  This  agreement  indicates  that  the  ideas  put  forth  in 
determining  the  distance  between  the  two  separations  are  along  the 
right  lines,  but  further  improvements  with  more  data  are  very  desira¬ 
ble.  It  should  be  noted  that  in  the  present  analysis,  the  possible  role 
of  another  geometrical  parameter  equivalent  to  the  wedge  angle  in 
Pearson’s  experiments  is  not  considered.  In  our  problem  the  above 
mentioned  parameter  would  be  equivalent  to  the  angle  of  boundary  layer 
separation.  Neglect  of  this  parameter  may  be  responsible  for  some  of 
the  scatter  in  the  observed  data  of  Fig.  15. 

A  straight  line  approximation  for  the  data  of  present  models  in 
Fig.  15  gives, 

<2> 

Since,  X  can  never  be  negative  it  should  be  taken  that 

^<Vn.c.’°  for  ¥  [1-'VN'c-]>'42 

and  for  values  of 
evaluate  X  . 


:p,)N.C.  ]  S  '42, 


equation  (2)  should  be  used  to 


6.  THE  COMPLETE  CORRELATION 


One  can  express  the  cavity  separation  position  as  the  sum  of 
the  position  of  the  boundary  layer  separation  preceding  it  and  the  dis¬ 
tance  between  the  two  or, 

<sr  *  /D)  =  (s„  T  q  /D)  +  (X*/D)  (3) 

using  the  correlations  obtained  earlier  in  equations  (1)  and  (2),  equation 
( 3)  may  be  written  as 


<8C.S./D»=  <sB.L.S./D>N.C.  -2'37  [I  +(tlN.C.]  * 

Ps 


['8B.L.S./D>N.C.-"'!pmin/D>N.C.;l+  <4> 

'VDkc.{130-J12lf[l-lcP,>K.C.1}' 

Therefore,  in  general 

<*C.S./D>  =  <8C.S./D>[R«D*  ^cps^N.C. '  We^*  ^ 

To  use  equation  (4),  the  knowledge  of  the  pressure  distribution  on  the 
headform  under  non-cavitating  conditions  is  a  requirement  and  the 
position  of  the  boundary  layer  separation  should  be  known  or  may  be 
estimated  from  the  method  outlined  in  Appendix  A.  The  method  of 
estimating  (6,)^  q  is  explained  in  the  same  Appendix.  It  should  be 
noted  that  equation  (4)  is  valid  only  for  -Hjr  Cl  -  (c  C  ^  <  *or 
values  of  this  term  greater  than  0.42,  the  third  term  vanishes. 


r~  f 
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7.  COMPUTATIONS  USING  THE  COMPLETE  CORRELATION 

Computations  were  carried  out  to  predict  the  position  of 
cavitation  separation  on  spheres  and  a  cylinder  for  various  Reynolds 
number  range  and  cavitation  numbers.  For  these  calculations  it  is 
convenient  to  convert  equation  (4)  in  terms  of  angular  positions,  then 

YC.S.=  (YB.L.S.*N.C.  “  2*37^1  +  (c7lN.C.^(YB.L.S.>N.C.  “  (Ycpmin*N.C.^ 


+  (w)  [4>N.C.  ]  {l30  -  312  M  [1-  <cp  )nx.  ]} 


(6) 


where  Y  represents  the  angular  position  from  the  stagnation  point.  The 
following  information  under  non-cavitating  conditions  was  obtained  with 
the  help  of  experimental  pressure  distributions  from  various  references, 

SPHERE 


a)  104<ReD<  105  -  From  Ref.  (51) 

<yb.l.s.>n.c.  =  82-5' 

(Ycpmin*N.C.  =  74° 

^p^N.C. =  ",4 


b)  Ren=  2  X  10  -  From  Ref.  (51) 

^B.L.S.^N.C.3  82,5 * 
^^min^N.C.572,50 

<cp>c-'46 


Below  the  first  critical 
Reynolds  number 


At  the  first  critical 
Reynolds  number 


m  -  suvaPVMftr  ** 
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c)  ReD>  3  x  103  -  From  Ref.  (51) 
(YB.L.S.*N.C.  =  109° 


(Ycpmin*N.C.=  90  ° 
^Cps^N.C,  * 


Above  the  first  critical 
Reynolds  number 


CYLINDER 

104<  R«d<  105  -  From  Ref.  (42) 


^YB.L.S.^N.C. =  80  * 
^Ycpmin^N.C. =  70* 

(cp>cLu875 


Below  the  first  critical 
Reynolds  number 


Momentum  thickness  at  the  boundary  layer  separation  under  non- 
cavitating  conditions  was  estimated  by  Thwaite's  method  and  is  shown 
in  Fig.  A.  5  of  Appendix  A. 

SAMPLE  CALCULATIONS  -  SPHERE 

U  =10  ft/sec.,  a=  .1,  D  =  l/2" 

Re^  —  4  x  104  from  Appendix  A  =»  (0g/D)jj  ^  —  14  x  10“4 

U=  2.05  X  10“5  -  .  T=.0047  ^ 

ft*  ft 

and  from  Ref.  (51 ),  as  mentioned  above 

^B.L.S.^N.C.”82,5*  *  ^Ycpmin^N.C.  ~74*  *  ^g^N.C.  ”  “  A 

*TT  Ci"(cp Jn.c.^  ~ *°6i  »  &  +  (c^N.C.^  =  *75, 

8  Ps 


Using  above  values  in  expression  (  6  )  one  gets 


Ycg  =82.5°-  [2.37  (,75)(82.5#  -  748)] 

+  {•—£■  X  14  X  10“4X  [130- 312  (.061)]} 

Yc  s  —  82.5®-  15°  +  17.8®  - 85.3°. 

Similar  calculations  for  sphere  and  cylinder  at  various  Reynolds 
number  and  cavitation  numbers  are  shown  in  Figs.  16,  17  and  18. 
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8.  COMPARISON  WITH  EXPERIMENTS 


Also  shown  in  Figs.  16,  17  and  18  are  experimentally  found 
positions  of  cavitation  separation  obtained  by  Brennen^’  Also 
marked  is  the  position  of  cavitation  separation  obtained  theoretically 

(53) 

by  using  the  smooth  separation  condition  and  computed  by  Brennen  . 

From  Fig.  16,  at  a=.l  the  calculated  values  follow  the  same 
trend  with  Re^  as  the  experimentally  observed  values  and  are  within 
7*.  From  Fig.  17,  at  o=  .25  the  agreement  is  much  better  and  the 
calculated  values  are  within  4°  of  that  found  experimentally  and  similar 
results  are  seen  for  o  =  0.5  in  the  same  Fig.  17.  The  important  point 
to  note  from  these  figures  is  thct  the  slope  of  the  calculated  curves 
and  experimentally  found  curves  is  in  excellent  agreement  for  all 
Reynolds  numbers  and  cavitation  numbers.  From  Fig.  17,  at  o=  0.25, 
the  discontinuity  in  the  observed  positions  just  before  and  after  the  first 
Reocrit  *or  a  8P^ere  *8  ai80  predicted  from  present  calculations. 

From  Fig.  18,  for  a  cylinder  the  agreement  is  very  good  at  an 

4  5 

Re^j  of  10  but  gradually  differs  by  about  10°  at  Re^  of  10  . 

Prediction  of  the  position  of  cavitation  separation  for  three 
different  fluids  which  possess  different  values  of  |i/T  is  shown  in  Fig. 
19.  The  value  of  the  coefficient  of  surface  tension,  T,  for  mercury  is 
about  8  times  that  of  water  and  for  benzene  it  is  one-half  of  that  for 
water.  It  may  be  concluded  from  Fig.  19  that  for  low  Reynolds  number 
range  (Re^  <8x10  ),  viscous  forces  dominate  in  determining  the 
position  of  cavitation  separation  on  a  sphere  and  this  is  consistent  with 
the  experimental  observations  of  Brennen^). 
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9.  CONCLUDING  REMARKS 


The  proposed  method  to  predict  the  position  of  cavitation 
separation  applies  only  to  bodies  which  possess  laminar  boundary 
layer  separation  under  fully  wetted  conditions. 

On  such  bodies  with  a  fully  developed  cavity,  cavitation  separa¬ 
tion  is  observed  to  be  preceded  by  a  laminar  boundary  layer  separation 
and  the  distance  between  the  two  separations  shows  a  strong  dependence 
on  Reynolds  number. 

The  present  method  is  based  on  data  acquired  on  only  two 
different  models;  further  acquisition  of  data  on  different  bodies  with 
present  experimental  techniques  would  strengthen  the  validity  of  and 
may  improve  the  proposed  empirical  method. 

The  predicted  values  of  the  position  of  cavitation  separation 
based  on  present  methods  on  a  sphere  for  Reynolds  number  range  of 
10^  to  10^  showed  good  agreement  with  the  experimentally  observed 
values  and  are  far  superior  to  that  predicted  by  the  smooth  separation 
condition  used  in  classical  free  streamline  analysis. 
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Fig.  1. 


Fig.  2. 


Fig.  3. 


Fig.  4. 


Fig.  5. 


Fig.  6. 


Fig.  7. 
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LIST  OF  FIGURE  CAPTIONS  -  PART  II 


Photographs  illustrating  the  striking  difference  in  the 
physical  appearance  of  cavitation  in  the  cases  of  nucleate 
cavitation  separation  and  viscous  laminar  cavitation  separation. 

Silhouette  schlieren  photograph  showing  existence  of  the 
laminar  boundary  layer  separation  upstream  of  cavitation 
separation  on  the  Swedish  headform  with  a  fully  developed 
cavity.  U=  25.6  ft/sec. ,  a  =  .39.  The  actual  distance  be¬ 
tween  the  two  separations  is  0.11  in. 

Variation  of  angular  position  of  cavitation  separation  with 
tunnel  velocity  and  cavitation  number  for  the  hemispherical 
nose. 

Variation  of  angular  position  of  boundary  layer  separation 
under  cavitating  conditions  with  tunnel  velocity  and  cavitation 
number  for  the  hemispherical  nose. 

Variation  of  the  streamwise  position  of  cavitation  separation 
on  the  Swedish  headform  with  tunnel  velocity  and  cavitation 
number. 

Variation  of  the  streamwise  position  of  the  boundary  layer 
separation  under  cavitating  conditions  on  the  Swedish  head- 
form  with  tunnel  velocity  and  cavitation  number. 

Variation  of  the  streamwise  distance  between  the  position  of 
boundary  layer  separation  and  the  cavity  separation  with 
tunnel  velocity  and  cavitation  number  for  the  hemispherical 
nose. 

Variation  of  the  streamwise  distance  between  the  position  of 
boundary  layer  separation  and  the  cavity  separation  with 
tunnel  velocity  and  cavitation  number  for  the  Swedish  headform. 


.  '  '"VW  «r>  , 


Fig.  9. 
Fig.  10. 
Fig.  11. 

Fig,  12. 
Fig.  13. 

Fig.  14. 
Fig.  15. 

Fig.  16. 

Fig.  17. 

Fig.  18. 

Fig.  19. 
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Data  of  Fig.  7  at  only  two  velocities  for  which  maximum 
amount  of  data  was  obtained. 

Data  of  Fig.  8  at  only  two  velocities  for  which  maximum 
amount  of  data  was  obtained. 

Variation  of  the  streamwise  distance  between  the  position  of 
boundary  layer  separation  and  the  cavity  separation  at  any 
cavitation  number  with  tunnel  velocity  for  both  the  test 
models. 

A  correlation  plot  to  determine  the  position  of  boundary  layer 
separation  under  cavitating  conditions. 

A  schematic  drawing  showingthe  apparatus  used  by  Pear  * 
son^*^  to  study  the  entrainment  of  an  air  cavity  in  the 
diverging  section  of  a  wedge. 

A  schematic  drawing  showing  the  physical  conditions  under 
which  cavitation  occurs  on  the  present  models. 

A  correlation  plot  to  determine  the  distance  from  the  position 
of  boundary  layer  separation  to  the  position  of  cavitation 
separation* 

Comparison  of  presently  calculated  position  of  cavitation 
separation  with  experimentally  observed  position  on  sphere 
for  10^<  10^  at  a=.l. 

Comparison  of  presently  calculated  position  of  cavitation 
separation  with  experimentally  observed  position  on  sphere 
for  10^<  10^  at  c  =  .25  and  .5. 

Comparison  of  presently  calculated  position  of  cavitation 
separation  with  experimentally  observe  i  position  on  cylinder 
for  10^<Rejj<10^  at  a  =.5. 

Comparison  of  calculated  position  of  cavitation  separation  on 
sphere  at  a=  .1  for  three  different  liquids  with  different 
physical  properties. 


Plate  1:  Cavitation  on  a  bi-  convex  hydrofoil  showing 
nucleate  cavitation  separation. 

U=  51  ft/sec  a  =.11  Ref.  6 


Plate  2:  Cavitation  on  Swedish  headform  showing 
viscous  laminar  cavitation  separation. 
(Courtesy  of  Dr.  A.  J.  Acosta) 

U=  40  ft/ sec  a  =.424  Ref.  9 


Fig.  1.  Photographs  illustrating  the  striking 
difference  in  the  physical  appearance  of  cavitation 
in  the  cases  of  nucleate  captation  separation  and 
viscous  laminar  cavitation  separation. 
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Fig. 3.  Variation  of  angular  position  of  cavitation  separation  wit] 
tunnel  velocity  and  cavitation  nxrrber  for  the  hemispherical  nose, 
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Fig.  10.  Data  of  Fig.  8  at  only  two 
velocities  for  which  maximum 
amount  of  data  was  obtained. 
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TUNNEL  VELOCITY  -U  IN  FEET  PER  SECOND 

Fig.  11.  Variation  of  the  streamvsise  distance  between  the  position  of 
boundary  layer  separation  and  tht  cavity  separation  at  any  cavitation 
number  with  tunnel  velocity  for  both  the  test  models. 
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Fig.  13.  A  schematic  drawing  showing  the  apparatus  used  by 
Pearson(SO)  to  study  the  entrainment  of  a  air  cavity  in  the  diverging 
section  of  a  wedge. 


r  separated 
streamline 


Fig.  14.  A  schematic  drawing  showing  the  physical  conditions  under 
which  cavitation  occurs  on  the  present  models. 
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REYNOLD'S  NUMBER  -ReD  =  UD/v 

Fig.  16.  Comparison  of  presently  calculated  position  of  cavitation 
separation  with  experimentally  observed  position  on  sphere  for 
10*<  Ren<  10®  at  cr=.l. 


3 


REYNOLDS  NUMBER- ReD  =  UDA 

.17.  Comparison  of  presently  calculated  position  of  cavitation 
aration  with  experimentally  observed  position  on  sphere  for 
<Rej)<10°at  a=.25and.5. 


CYLINDER -0=1/4",  <r  =  0.5 


Fig.  18.  Comparison  of  presently  calculated  position  of  cavitation 
separation  with  experimentally  observed  position  on  cylinder  for 
10*<Ren<105  at  ct=  .5. 


"I 

OT  10 

'  I 


ro  * 


s  a 

7  2 


w  NJ 
H  2 


o 

<7> 

s 

>- H 
t  uj 
oo 

8| 

u.  u. 

o  o 

I-  H 

z  z 

UJ  UJ 

o  o 

u.  il 
u.  u. 
uj  uj 
O  O 
o  o 

h  n 


w>,e  o 
JH  P* 
n  io  p. 


SS3M93Q  Nl  S'0/-3"I9NV  N0liV«Vd3S  NOIlVilAVO 


REYNOLDS  NUMBER  -  Re0  x  I0~H 

;.19.  Comparison  of  calculated  position  of  cavitation  separation  on 
tere  at  CT  =.  1  for  three  different  liquids  with  different  physical 
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APPENDIX  A 

LAMINAR  BOUNDARY  LAYER  GROWTH  CALCULATIONS 


The  laminar  boundary  layer  growth  calculations  were  made  by 
using  Thwaites  approximate  method  which  is  detailed  in  Ref.(54)  for 
two  dimensional  bodies  and  is  detailed  in  Ref.  (55)  for  axisymmetric 
bodies.  Computations  were  made  on  a  hemispherical  nose,  Swedish 
headform,  1.5  calibre  ogive,  a  sphere  and  a  cylinder. 

Necessary  information  to  use  the  Thwaites  method  is  the  knowl¬ 
edge  of  the  pressure  distribution  around  the  headform.  Presently  used 
experimental  pressure  distributions  obtained  from  different  sources 
are  shown  in  Figs.  A.  1  and  A.2.  It  is  assumed  that  the  gross  pressure 
distribution  around  the  hemispherical  nose,  the  Swedish  headform  and 
the  1.5  calibre  ogive  does  not  change  with  Re^  for  Re^  >  2  X  10*\  This 
assumption  may  be  justified  from  the  observations  of  Part  I,  that  the 
hemispherical  nose  and  the  Swedish  headform  possess  very  small 
laminar  separated  regions  for  ReD  >  2  X  10**.  Rouse^3^  has  observed 

that  the  gross  pressure  distribution  around  the  hemispherical  nose  does 

5 

not  change  appreciably  with  Re^  for  Re^  >2X10. 

In  Thwaites  method  a  parameter  X  is  defined  as  follows, 


where  0  is  the  momentum  thickness,  U#  is  the  local  velocity  gradient 
at  the  edge  of  the  boundary  layer  and  v  is  the  kinematic  viscosity.  Then, 
from  Thwaites  approximate  method, 


R 
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„AS  <  f 
'  2  ?.**' 


2„5  , 
rAU  dx 
0  e 


(A.  2) 


where  s  is  the  streamwise  coordinate  along  tL  model  surface,  is 
the  local  radius  as  shown  in  Fig.  A.  3  and  U  is  the  local  velocity  at 
the  edge  of  the  boundary  layer.  Non-dimensionalizing  velocities  with 
which  is  the  velocity  far  away  from  the  body  and  distances  with 
diameter  of  the  body,  D,  in  equation  (A. 2)  one  obtains 


=  J£[.{dc  M.,!,}.  1_], 


i8>^)  5/2 

I  r0tl“cn-1  dx-  (A*3) 


U_ 


Also  used  in  equation  (A.  2)  is  the  relationship  =  i/l^t  where  c  is 
the  local  pressure  coefficient  defined  by  (P-P  1/1/2  pT  .  P  is  the 

00  Xj  00  00 

pressure  far  away  from  the  body  where  velocity  is  and  p^  is  the 
liquid  density. 

The  integration  of  equation  (A.  3)  was  carried  out  by  trapezoidal 
rule  with  the  help  of  known  experimental  pressure  distribution.  Thus, 
once  X  is  determined  from  equation  (A.  3)  the  shape  parameter,  H, 
which  is  the  ratio  of  displacement  thickness  to  momentum  thickness, 

jjC 

6/0,  may  be  found  from  page  306  of  Ref.  (55).  The  modified  values  of 
Curie  and  Skan  were  used. 

Knowing  X,  one  may  determine  0  and  H  which  provides  the 
knowledge  of  6  .  Laminar  boundary  layer  separation  is  said  to  occur 


when  X  =  -0.09  with  corresponding  H=  3.55. 
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Results:  Geometrical  details  of  a  hemispherical  nose  and  Swedish 

headform  are  shown  in  Figs.  9  and  10  of  Part  I  and  that  of  1.5  calibre 
ogive  are  shown  in  Fig.  A.4.  Comparison  of  the  predicted  position  of 
the  laminar  boundary  layer  separation  by  Thwaites  approximate  method 
for  several  cases  with  that  observed  experimentally  is  shown  in  Table 
A.  1.  The  predictions  are  seen  to  be  within  7  percent  of  experimental 
observations.  Variation  of  the  momentum  thickness  at  the  predicted 
position  of  separation,  P  ,  with  Reynolds  number  for  different  bodies 
is  shown  in  Fig.  A.  5. 


TABLE  A.l  Comparison  of  Calculated  Positions  of  Boundary  Layer  Separation  by  Thwaites’  Method  with 
Experimentally  Observed  Positions  for  Various  Bodies. 
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Fig.  A.  1. 

Fig.  A.  2. 

Fig.  A.  3. 
Fig.  A.4. 

Fig.  A.  5. 


LIST  OF  FIGURE  CAPTIONS  -  APPENDIX  A 


Pressure  distribution  on  the  hemispherical  nose,  the  Swedish 
headform  and  the  1.5  calibre  ogive.  Taken  from  Refs.  (4)  and 
(43). 

Pressure  distribution  on  a  sphere  and  a  cylinder  in  various 
Reynolds  number  ranges.  Taken  from  Refs.  (42)  and  (51). 

Coordinate  system  used  for  boundary  layer  growth  calculations. 

Diagram  showing  the  geometrical  details  of  the  1.5  calibre 
ogive. 

Variation  of  calculated  value  of  momentum  thickness  at 
separation  by  Thwaites*  method  with  Reynolds  number  for 
various  models. 
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Fig.A.2.  Pressure  distribution  on  a  sphere  and  a  cylinder  in 
various  Reynolds  number  ranges.  Taken  from  Refs.  (42)  and  (51). 


Fig,  A. 4.  Diagram  showing  thegeometrical  details  of  the  1.5 
calibre  ogive. 
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Fig.  A. 5.  Variation  of  calculated  value  of  momentum  thickness  at 
separation  by  Thwaites  method  with  Reynolds  number  for  various 
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APPENDIX  B 

APPROXIMATE  COMPUTATIONS  OF  SPATIAL  GROWTH  OF 
BOUNDARY  LAYER  DISTURBANCES  TO  PREDICT  TRANSITION 

In  Appendix  B  an  approximate  method  of  predicting  transition  of 
the  laminar  boundary  layer  to  turbulent  boundary  layer  using  the  spatial 
stability  charts  for  the  Falkner-Skan  boundary  layer  profiles  is  pre¬ 
sented.  The  stability  charts  have  been  computed  by  Wazzan,  Okamura, 
and  Smith  and  are  presented  in  Refs.  (46)  and  (41).  An  approximate 
method  of  obtaining  necessary  boundary  layer  growth  information  for 
the  use  of  stability  charts  is  outlined  in  Appendix  A.  Transition  of 

boundary  layer  is  said  to  occur  when  the  amplitude  of  initial  disturbance 

7 

is  magnified  by  a  factor  of  e  ,  which  is  based  on  the  semi -empirical 

method  of  Smith  (41).  For  axisymmetric  bodies  the  calculated  amplifi- 

7  9 

cation  factor  by  Smith  was  found  to  be  about  e  toe  at  the  experimental 

position  of  transition.  This  method  has  been  applied  for  hemispherical 

nose,  Swedish  headform  and  1.5  calibre  ogive.  The  value  of  Reynolds 

number  at  which  transition  occurs  at  the  position  of  laminar  boundary 

layer  separation  is  said  to  be  a  critical  value  and  is  designated  by 

]%eDcrit’  Roshko  (42)  refers  to  this  Reynolds  number  as  second 

critical  Reynolds  number.  It  should  be  noted  that  for  Rert<  Re_ 

u  u  crit 

transition  to  turbulence  occurs  on  the  separated  free  shear  layer  and 
thus  cannot  be  predicted  by  the  present  method. 

In  the  method  of  Ref.  (46)  the  stream  function  is  assumed  to  be 


of  the  form 
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♦0(y)+*(*.y.t)  (B.i) 

where  represents  the  mean  flow  and  t  a  small  disturbance  commonly 
known  as  Tollmein-Schlichting  disturbance  of  the  form 

*  =  cp(y)ei<X(x‘ct)  (B.  2) 

In  the  above  two  equations  the  distances  are  normalized  by  boundary 

layer  thickness  5,  velocity  by  U  ,  velocity  at  the  edge  of  the  boundary 

layer  and  t  by  5/U  .  Substitution  of  these  equations  in  continuity  and 

Navier -Stokes  equations  of  motion  and  linearizing  them  in  cp  produces 

the  Orr-Sommerfield  equation.  It  should  be  ncted  that  a  in  general  is 

complex  and  imaginary  part  ou  is  the  disturbance  amplification  rate  and 

the  real  part  a  the  wave  number,  and  c  is  complex  of  which  the  real 
r 

part  cr  is  the  phase  speed.  Since  only  spatial  amplifications  are 
considered  the  product  ttc  is  taken  to  be  real  quantity.  The  Orr- 
Sommerfield  equation  with  proper  boundary  conditions  is  solved  by 
numerical  means.  The  details  should  be  referred  to  in  Ref.  (46).  The 
Falkner-Skan  velocity  profiles  are  obtained  by  numerical  integration  of 
the  boundary  layer  equation 

F*(T))  =  -FF#  +  p(F'Z-l)  (B.3) 

where  ri  =  y/6,  F^y)  =  U(y)/U  and  prime  denotes  differentiation  with  r\. 

It  should  be  remembered  that  (3  >0  represents  accelerating  flow,  (3  =  0 
represents  flat  plate  solution  and  (3<0  represents  deaccelerating  flow. 

(3  =  -.1988  represents  separating  velocity  profile.  Therefore  solution  of 
Eq.  (B.  3)  provides  a  value  for  the  shape  parameter  H=  6*/0  as  function 
of  (3.  Tabulation  of  values  of  H  for  different  (3,s  is  given  on  page  20  of 
Ref.  (46)  and  are  shown  in  Fig.  B.I . 
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The  stability  charts  of  Ref.  (46)  are  presented  in  the  form  of 

curves  of  constant  a.  6  /R  .  on  the  diagrams  of  dimensionless  frequency 

x  g* 

vs.  Reynolds  number  {u)  ,  R  *).  This  is  done  for  many  values  of  p. 

r 

These  charts  represent  the  eigenvalue  solutions  of  Orr-Sommerfield 
equation. 

If  8.  is  taken  to  be  the  neutral  point  or  position  from  where  the 
disturbances  begin  to  grow,  then  amplification  at  any  point  s  downstream 
of  s.  in  terms  of  tabulated  quantities  is  given  by 

Ps/Q  x  U 

£nA(s/D)  =  -J  (ai6*/Re  ^  .  ReD .  d(s'/D)  (B.4) 


The  above  equation  is  for  fixed  U)  =  0) *v  /U?  . 

re 

Method  of  Present  Computations:  From  Thwaites'  method  local  value  of  Rg 

H  as  function  of  s/D  is  calculated,  corresponding  to  this  H  from  Fig.  B.l 

a  p  is  ob<  lined.  At  this  p  corresponding  chart  in  Ref.  (46)  gives  the  value 

Ue 

of  (a.  6*/R  *)  for  different  value  of  U)  .  Local  value  ofw—  is  obtained 
1  6  r  Uoo 

from  experimental  pressure  distribution.  Many  times  extrapolation  of 

charts  from  Ref.  (46)  had  to  be  carried  out  by  replotting  the  results. 

*  Ue 

Therefore  the  product  of  (ctjft  /Re  „,)  •  —  •  Re^  can  be  plotted  as  function 

ft  co 

of  s/D.  This  has  to  be  done  for  different  values  of  U'r>  and  among  these 
u>r  one  particular  value  will  give  the  maximum  amplification  at  the 
desired  position.  Thus  fn(Amax)  may  be  calculated  as  function  of  posi¬ 
tion  along  the  model  surface  by  simple  integration.  Such  computations 
for  hemispherical  nose  are  shown  in  Fig.  B.  2;  it  may  be  seen  that  even 
at  a  Reynolds  number  of  2.2  *10^  transition  has  not  occurred  prior  to 


-145- 


predicted  boundary  layer  separation.  The  critical  Reynolds  for  hemi¬ 
spherical  nose  was  estimated  to  be  about  5  XIO^.  It  was  seen  earlier 

that  during  present  experiments  schlieren  method  showed  a  laminar 

5 

separation  bubble  to  exist  at  a  Reynolds  number  of  9  X  10  confirming  the 
predictions  afforded  in  Fig.  B.2.  Computation  of  maximum  amplifica¬ 
tion  rates  for  the  Swedish  headform  are  shown  in  Fig.  B.3.  The  critical 
Reynolds  number  for  Swedish  headform  is  found  to  be  around  1.5X10^ 
and  from  schlieren  experiments  a  laminar  separation  bubble  existed  on 

5 

the  Swedish  headform  at  a  Reynolds  number  of  8.78X10  .  For  Rep> 

1.5X10^  transition  is  predicted  prior  to  the  predicted  and  observed 

position  of  laminar  boundary  layer  separation  as  shown  in  Fig.  B.3. 

From  Fig.  B.4,  the  estimated  critical  Reynolds  number  for  1.5  calibre 
c 

ogive  is  6.28X10  .  Also  shown  in  Fig.  B.4  are  predicted  position  of 

transition  for  various  Rep  greater  than  6.28X10  .  Once  knowing  the 

position  of  transition,  one  may  estimate  the  magnitude  of  the  pressure 

coefficient  at  transition,  c_  ,  from  the  experimentally  measured  pres- 

ptr 

sure  distribution.  Using  the  results  of  Fig.  B.4  in  conjunction  with 


pressure  distribution  for  1.5  calibre  ogive  shown  in  Fig.  A.  1  of 


Appendix  A,  variation  of  -Cp^  with  Rep  is  shown  in  Fig.  B.5.  As 
expected  the  position  of  transition  moves  forward  with  increas  in  Rep 


and  hence,  -Cp^  may  be  expected  to  increase  with  increase  in  Rep  as 
observed  in  Fig.  B.5. 


A  few  concluding  remarks  should  be  made  concerning  the  method 
outlined  above.  During  the  computations  described  earlier,  an  assump¬ 
tion  has  been  made  that  the  gross  pressure  distribution  around  the 
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5 

headform  does  not  change  with  Re^  for  Re^>  2  X10  .  This  assumption 
has  been  justified  in  Appendix  A.  It  should  be  noted  that  for  axi- 
symmetric  bodies  the  difference  on  amplification  rates  from  that  of  a 
two  dimensional  body  comes  about  only  through  the  change  of  boundary 
layer  growth  properties. 


A, 
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Fig.  B.l. 
Fig.  B.2. 

Fig.  B.3. 

Fig.  B.4. 

Fig.  B.5. 


LIST  OF  FIGURE  CAPTIONS  -  APPENDIX  B 


Relationship  between  Hartee  (3-dimensionless  pressure 
coefficimt  and  H- shape  parameter.  Taken  from  Ref.  (46). 

Calculated  maximum  amplification  factors  along  streamwise 
position  of  the  hemispherical  nose  with  increasing  Reynolds 
number. 

Calculated  maximum  amplification  factors  along  streamwise 
position  of  the  Swedish  headform  with  increasing  Reynolds 
number. 

Calculated  maximum  amplification  factors  along  streamwise 
position  of  the  1.5  calibre  ogive  with  increasing  Reynolds 
number. 

Variation  of  the  calculated  pressure  coefficient  at  the  pre¬ 
dicted  position  of  transition  with  Reynolds  number  for 
Rejj>Rejjcrjt.  The  results  are  for  the  1.5  calibre  ogive. 
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Fig.  B.2.  Calculated  maximum  amplification  factors  along  streamwise 
position  of  the  hemispherical  nosewith  increasing  Reynolds  number. 
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Fig.  B.3.  Calculated  maximum  amplification  factors  along  streamwise 
position  of  the  Swedish  headformwith  increasing  Reynolds  number. 
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Pig.  B.4.  Calculated  maximum  amplification  factors  along  streamwise 
position  of  the  1.5  calibre  ogive  with  increasing  Reynolds  number. 
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Fig.B.5.  Variation  of  the  calculated  pressure  coefficient  at  the 
predicted  position  of  transition  with  Reynolds  number  for 
ReD>ReD  crJt.  The  results  are  for  the  1.5  calibre  ogive. 
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